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ABSTRACT 

Black hole masses predicted from the M, — a relationship conflict with those pre- 
dicted from the M, — L relationship for the most luminous galaxies, such as brightest 
cluster galaxies (BCGs). This is because stellar velocity dispersion, a, increases only 
weakly with luminosity for BCGs and other giant ellipticals. The M, — L relationship 
predicts that the most luminous BCGs may harbor black holes with M, approaching 
10^^ Mq, while the M, — a relationship always predicts M, < 3 x 10^ Mq. Lacking direct 
determination of M, in a sample of the most luminous galaxies, we advance arguments 
that the M, — L relationship is a plausible or even preferred description for BCGs and 
other galaxies of similar luminosity. Under the hypothesis that cores in central stellar 
density are formed by binary black holes, the inner-core cusp radius, r-y, may be an 
independent witness of M,. Using central structural parameters derived from a large 
sample of early-type galaxies observed by HST, we argue that L is superior to a as an 
indicator of r-y in luminous galaxies. Further, the observed r-y — M, relationship for 11 
core galaxies with measured M, appears to be consistent with the M, — L relationship 
for BCGs. BCGs have large cores appropriate for their large luminosities that may be 
difficult to generate with the more modest black hole masses inferred from the M, — a 
relationship. M, ~ L may be expected to hold for BCGs, if they were formed in dissipa- 
tionless mergers, which should preserve ratio of black hole to stellar mass. This picture 
appears to be consistent with the slow increase in a with L and the more rapid increase 
in effective radii, Re, with L seen in BCGs as compared to less luminous galaxies. If 
BCGs have large BHs commensurate with their high luminosities, then the local black 
hole mass function for M, > 3 x W^Mq may be nearly an order of magnitude richer 
than what would be inferred from the M, — a relationship. The volume density of 
the most luminous QSOs at earlier epochs may favor the predictions from the M, — L 
relationship. 

Subject headings: galaxies: nuclei — galaxies: structure — black hole physics 
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The Most Luminous Galaxies 



The Most Massive Black Holes 



Nearly every elliptical galaxy and spiral bulge has a black hole at its center (jMagorrian et al 



19981 ). The masses of the black holes, M,, are r elated to the F-band luminosity, L, and average stel- 
lar velocity dispersion, a, of their host galaxies (lDressleJl989l:lKormendvlll993l :l iKormendv &: Richst 



1995 



Magorrian et al.lll998l : lFerrarese &: Merritt 



Haring fc Ri3ll2nn4l ). The M, 



2001 



I 



Gebhardt et al 



2000a 



Tremaine et al 



one 



20021 : 



prediction of black hole masses 
galaxy parameters. 



a and M, — L relationships are powerful tools as they allow the 
— which are difficult to measure directly — from readily available 



The black hole population in the most massive galaxies has yet to be assayed, however, which 
means that estimates of M, in these objects are based on extrapolations of relationships defined 
by smaller galaxies. The current reco rd for largest black hole mass measured directly is M, ~ 
3 X 10^ Mq in M87 ( Harms et al.l 11994 ). yet M87 is only the second-ranked galaxy in a cluster of 
modest richness . Brightest cluster gala xies (BCGs) in nearby Abell clusters are typically ~ 3x 
more luminous (jPostman &: Laueiill995l ) and may host proportionately more massive BHs. Testing 
this hypothesis through measurements of stellar dynamics requires both high sensitivity and high 
spatial-resolution, given the low central surface brightnesses and relatively large distances of BCGs. 
Such observations were not possible with the Hubble Space Telescope (HST) even before the failure 
of the Space Telescope Imaging Spectrograph; they are only now becoming feasible with the advent 
of adaptive optics spectroscopy on 10m class telescopes. 

A number of arguments suggest that black holes w it h M« > 3 x W^Mq do exist, even if 
this conclusion is not universal (e.g. iMcLure et al.l 12004 ). iNetzed (|2003l ) argues that some QSOs 
have M, > lO^^M© b ased on an emp i rical relatio nship between M,, broad-line width and nuclear 
luminosity for AGN. iBechtold et al.l (|2003l ) and IVestergaardI (|2004l ) also argue that some QSOs 
have black holes approaching this mass. Of particular relevance for BCGs is the hypothesis that 
cluster cooling flows a re inhibited by AGN heating from the central galaxy (jBinnev &: Taboiill995l : 
Churazov et al.l |2002| ) . Recent Chandra observati ons support a picture in which episodic AGN 
outbursts in BCGs heat the intra-cluster medium (jVoit &: Donahuell2005l): the energe tics required 
to terminate cooling flows imply M, > W^^Mq for many clusters ( Fabian et al.ll2002l ). 



Arguments for such massive black holes appear to be in conflict, however, with the expectations 
from the M, — a relationship applied to the local galaxy velocity-dispersion distribution function. 



^Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. These observations are associated with GO and GTO proposals ^ 5236, 5446, 5454, 5512, 
5943, 5990, 5999, 6099, 6386, 6554, 6587, 6633, 7468, 8683, and 9107. 

^The National Optical Astronomy Observatory is operated by AURA, Inc., under cooperative agreement with the 
National Science Foundation. 
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Tremaine et alJ (12002) find 

log(M./M0) = (4.02 ± 0.32) log(f7/200 km s~^) + 8.19 ± 0.06, 
for Hq = 70 km s^^Mpc^^ (which we will use throughout this paper). The lSheth et al.l ( 20031 ) local 



(1) 



velocity dispersion function shows a strong cut-off at o" ~ 400 km s~^ ^ which implies that ga laxies 
harboring black holes with M, > 3 x 10^ Mq would be extremely rare. iBernardi et al.l ( 2006al ) have 
identified a handful of galaxies with a > 400 km s^^, but their results do not alter this conclusion. 

Extrapolation of the M, — a relationship to galaxies more massive than M87 assumes that 
a (and not galaxy mass) is the fundar nental p a rame ter for determining M,. The uncertainty in 
such an extrapolation is underscored by I Wyithd (120061 ). who argues that the M, — a relationship is 
curved rather than linear in log-log space, in the sense that, at the high-c end, the "log-quadratic" 
relationship predicts higher M, than does equation ([T]). The Wyithe M, — a relationship, implies 
that the space density of black holes with M, > 5 x 10^ Mq may be substantially higher than that 
implied by equation ([1]) (although the exact difference is highly sensitive to both the details of the 
velocity dispersion distribution function, and the assumed level of cosmic scatter in the M, — a 
relationship) . 

In this paper we point out that the M, — L relationship applied to the most lu minous galaxies 



predic ts M, values tha,t are significantly larger than those predicted by either the iTremaine et al 



( 2002) or IWvithj l200 d ) M, — a relationships. This difference arises because BCGs do not follow the 
Faber &: JacksonI (119761 ) relationship between L and a. The relationship be tween L and a "plateaus" 



at large L in the sense that BCGs have relatively low a for their high L (lOegerle &: Hoessellll991 



see also iBovlan-Kolchin et al.ll2006l . who have seen this effect in simulations.) 



Resolution of which of the M, — L or M, — a relationships is most representative of the black 
hole population in the most massive galaxies will only be possible when black hole masses can 
be measured in such galaxies. In advance of such work, however, we can advance a number of 
arguments that suggest that the M, — L is a plausible and perhaps even preferred description for 
such systems. 

The first set of arguments are based on the central structure of BCGs and other luminous 



ellipt ical galaxies that have cores in their central brightness profiles (jLauer et al.lll995l : iLaine et al 
20021 ). A core is evident as a radius at which the steep envelope of the galaxy "breaks" and 



transitions to an inner cusp with a shallow slope in logarithmic coordinates. The favored theory 
for core formation posits that cores are formed wh en stars are ejected fr o m the galaxy's center 



by the decay of a. binary BH created in a merger (IBeeelman et al 



Faber et al.l 119971 : iQuinlan h HernquistI 119971 : iMilosavlievic Merrittll200ll ). The size of the core 



198nl: lEbisuzaki et al.1 [l991 



then reflects the total mass ejected, which should be a function of M,. The size of the core may 
thus be an independent witness of M,. In BCGs and other galaxies of similar luminosity, galaxy 
luminosity is more closely related to the physical scale of the cores than o", and the observed core 
size M, relationship for galaxies with cores and directly measured black hole masses appears to be 
consistent with the M, — L relationship. 
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A second set of arguments come from considering the formation of BCGs. If BCGs are formed 
in "dry" mergers, then the ratio of black hole to stellar mass should be preserved over mergers, 
leading to the observed M, — L relationship. In contrast, a may change little over such mergers, 
and no longer track black hole mass as well it does for the less luminous galaxies from which the 
M, — a has been determined. 

Lastly, we consider the relative predictions of the M, — L and M, — a relationships for the 
volume mass distribution function of black holes and which we compare to the predictions from 
QSO luminosity functions. A decisive discrimination between the two relationship s is not possible 



witho ut a better understanding of the cosmic scatter in both relationships, but the lTremaine et al 



(|2002l ) version of the M, — L relationship probably predicts too few high mass black holes to support 



the QSO luminosity function. 



2. A Large Sample of Early Type Galaxies With Central Structure Characterized 

by HST 

We start by comparing the two predictions M,(o") (Af, predicted from the M, — a relationship) 
and M,{L) (M, predicted from the M, — L relationship) for a sample of 219 gal axies for which we 
have central structural parameters derived from iJS'T imagery (jLauer et al.ll2007al ) . We then present 
the separate relationships between core structure versus a and L. This leads in turn to two separate 
predictions for how core size should be related to M,, which can be compared to the observed 
relationship between core size and M, for 11 core galaxies that have direct M, determinations. 

The galaxy sampl e combines several different HST imaging programs that all used the Nuker- 
law parameterization (jLauer et al.l Il995l ) to characterize the cen tral starlight distributions. The 



properties and definition of this sample are p resented in d e tail in 
we combine surface pho t ometry presented in lLauer et al.l (119951). 



iauer et al 



Faber et al.l (Il997l'l. 



(l2007al ). but briefly, 



Laine et al 



(120021 ). iRest et al.l (1200 il l. iRavindranath et al.l (120011 ). lOuillen et al.l(|200m . and lLauer et al.l (|2005l ). 

This diverse source material has been transferred to a common photometric system (y— band) and 
a common dis tance scale, adoptin g Hq = 70 km s~^ Mpc~^. The primary source of dis tances is the 



SBF s urvey of lTonry et al.l (|200ll ). but when possible we use the group mern ber ships in lFaber et al 



(jl989l ) and average SBF distances over the group. As the lTonry et al.l (j200ll ) SBF scale is consistent 
with Hq = 74, we sc ale up their S BF dis tances by 6%. The treatment of galaxies not in the SBF 
survey is detailed in iLauer et al.l (2007a). The sample is listed in Table [H It comprises 120 core 
galaxies, 87 power-law galaxies, and 12 intermediate galaxies. 

The most important Nuker-law parameter for the present analysis is the break radius, r^, which 
is used to calculate the cusp radius, r^, which in turn is used to represent the physical scale of the 
core (this parameter is discussed in detail in §4. II and Appendix IC]) . The average error in is 30%, 
based on comparison of Nuker parameters to non-parametric estimates of the same parameters. 



Central velocity dispersions are provided by the "Hyperleda" augmentation of the lPrugniel Simien 
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(|l996l ) compendium of published velocity dispersions; no values were available for 30 of the total of 
219 galaxies. We adopt a 10% typical error in a. The M, — a relationship as initially presented by 



Gebhardt et alj (j2000al ) used the average luminosity-weighted velocity dispersion measured in a slit 



along the major axis interior to the effective radi us. Velocity disp e rsion p rofiles are unfortunately 
not available for the bulk of the galaxies; however, iGebhardt et al.l ( 2000al ) showed that the central 
values are likely to be within 5% of the radial averages. 



2.1. Galaxy Luminosities 



1991 



The sources of the present galaxy luminosities are discussed in detail in lLauer et al.l (l2007al). 
Most of the magnitudes are derived from Vt or Bt values drawn from the RC3 (jde Vaucouleurs et al. 
). Bulge luminosities are given for SO and spiral galaxies base d on bulge/disk decompositions 



in the literature. Absolute luminosities were calculated using the ISchlegel et al.l (jl998l ) Galactic 
extinction values; we assume a typical My error of 10%. 

The accuracy of the BCG luminosities is of special concern as we will argue that they imply 
higher M, than would be inferred from the a values for the same galaxies. The present BCG 
luminosities are based o n fitting r^/^ laws to the inner portions (r < 50 kpc) of the R-band 
Postman &: Laueii (119951) brightness profiles, limiting the fits to radii that are well matched by 



this function. iGraham et al.l (|l996l ) show that BCG brightness profiles are better described by 
Sersic profiles with S ersic n > 4, which is also true of giant elliptical galaxies in general (e.q. 



Ferrarese et al. 



2006 



Kormendv et al.l 120071 ). However, BCGs with Sersic n > 4, typically have 



extremely large effective radii that are factors of several larger than the actual radial limit of the 
surface photometry; this in turn implies unrealistically large total luminosities. The r^/^ laws 
give a conservative lower limit for BCG total luminosities. Even so, the derived luminosities are 
systematically much larger than those provided by the Sloan Digital Sky Survey (SDSS). We resolve 
this issue in Appendix |A] with a demonstration that the SDSS BCG luminosities are strongly biased 
to low values by excessive sk y subtraction. The NIR ap parent magnitudes provided by the 2MASS 
Extended So urce Catalogue (IJarrett et al.ll2000l . 120031 ) have also been used to provide BCG total 
luminosities ([Batcheldor et al.ll2006l ): however in Appendix IB] we show that the 2MASS apparent 
magnitudes are also likely to be underestimates. 

A separate issue raised by a number of our colleagues is that BCG luminosities may need 
to be "corrected" for intracluster light (ICL). One such treatment of ICL assumes that the BCG 
is coincident with the center of the cluster p otential, and that th e composite BCG+ICL can be 
modeled as two superimposed r^/^ laws (cf. IConzalez et al.ll2005l ). The ICL component is then 
subtracted to yield the "true" BCG luminosity. A key feature of such models is that the ICL profile 
is assumed to continue to rise in brightness at radii well interior to where it dominates, thus implying 
a substantial contribution at even small radii. There is little physical justification for a correction 
of this form, however. As noted above, giant elliptical galaxies in general (not just BCGs) have 
Sersic n > 4. Further, the presumption that BCGs sit exactly at the center of their clusters is an 
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idealization that is actually realized in only a small fraction of systems. iPostman &: Laueii (|l995l ) 
show that BCGs are typically d isplaced frorn the geometric cluster center by ~ 90 kpc in projection 
and ~ 260 km in velocity. iPatel et al.l ( 20061 ) showed that BCGs ar e typically displaced frorn 
the centroid of cluster X-ray emission by 129 kpc, consistent with the iPostman &: Laueii (|l995l ) 
analysis. Lastly, the presumption that the ICL follows an r^^^ law into small radii is not uniquely 
demanded, and is probably inconsistent with the large velocity dispersion of stars truly not bound 
to t he BCG. Again, BCG s are well described over a large radial range by Sersic laws; in no case in 
the lGraham et al.l (119961 ) sample are there any profiles that have a distinct feature that objectively 
supports a two component model. This is not to say that ICL is not present, but the surface 
brightness at which it dominates even in the two component models are well outside the radii at 
which we me a sure t he r^/^ laws used to estimate total luminosity (typically less than 50 kpc). The 
Zibetti et al.l (|2005l ) models of ICL show that it begins to dominate the BCGs at r ~ 80 kpc from 
the BCG centers, corresponding to fir ~ 26. We conclude that a strong correction to our BCG 
luminosities for ICL is poorly justified. 



3. A Contradiction Between the M, — a and M, — L Relationships 



The M, — L rel ationship emerged from the first attempts to relate black hole ma ss to properties 
of the host galaxy (iDressleii Il989l : iKormendyl Il993l : iKormendy &: Richstond Il995l ) . Much of the 
recent work on this problem , however, has foc used on the M, — a relationship due to its apparent 
smaller scatter (although see lNovak et al.ll2006l on the significance of this), as well as arguments that 
a, rather than gal axy luminosity is the more fundamental parameter that determines how galaxies 



were formed (e.g.. lWvithe &: Loebll2005l ). While L and a are related by the lFaber &: JacksonI (|l976l ) 



relationship, since the disco very that galaxies lie on a "fundamental-plane " determined by L, a, 
and the effective radius, Re, (IDjorgovski Daviall987l : iDressler et al.lll987l ). we know that neither 
L nor 0" alone is sufficient to codify the full range of galaxy properties. The M, — L relationship 
thus may contain information that is not a trivial projection of the M, — a relationship. 

The relationsh ip between M, and L is shown in Figure [TJ Most of the galaxies shown are 
those presented in iTremaine et al. I (|i002)|f| transformed to Hq = 70 km s ^ Mpc ^. Due to the 
large scatter of the data points in Figure [H estimating a mean M, — L relationship is likely to be 
se nsitive to th e fittin g algorithm. We have elected to use the "symmetric" least-squares algorithm 
of iPress et al.l (|l992l ) throughout this analysis. This technique allows for errors in both variables 
being fitted, and finds the best slope and intercept parameters without assigning either parameter 
as the independent or dependent variable. As a way of bracketing uncertainties in the mean M, — L 



^We aug ment the Tremaine et al. (2002) sa mple with M, dete rminations in NGC 1399 (Houghton et al. 20061 ) ■ 
NGC 3031 jBower et al.l boOOh . NGC 3998 jBower et al.l boool ). NGC 4374 JSower et all Il998l ). NGC 4486B 
( Kormendv et al .I ll997l}. NGC 4945 dGreenhill et al.lll997l ). NGC 5128 (|Marconi et al.ll200ll ). NGC 7332 jNelson et al " 
2000l ). and Cygnus A (jTadhunter et al.ll2003l ). 
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relationship, we performed one fit using all the data points, but for a second fit we used only galaxies 
with My < —19, because they appear to have less scatter. The fit to all data points gives 

log(M./MQ) = (1.40 ± 0.17)(-My - 21)/2.5 + 8.41 ± 0.11, (2) 

which is shown as the dashed line in Figure [H Just fitting galaxies with My < — 19 gives 



log(M./M0) = (1.70 ± 0.22)(-My - 21)/2.5 + 8.22 ± 0.08, 



(3) 



which is shown as the dotted line in Figure [TJ Both relationships agree well for — 23 < My < 
—19; their differences in slope cause them to diverg e slightly when extra polated to more luminous 
galaxies. Both relationships also agree well with the lHaring &: Rixl (120041 ) relationship between M, 
and galaxy mass tr ansformed back to luminosity, which we consider as a third M, — L relationship. 
Novak et al.l (|2006l ) found that the M,-mass relationship was not significantly less tight than the 
M, — a relationship, given the errors of the various samples. If so, then the reduced scatter in 
the M,-mass relationship means that it should serve well as a relationship between M, and L; we 
tr ansform it by adopting M/Ly 6 x lQ~^-^^'^^^^~^'^'^^ Mq/ Lq, based on the M/L estimates given 



m 



Gebhardt et al.l (|2003l ): this gives 



log(M,/MQ) = (1.38 ± 0.07)(-My - 22)/2.5 + 8.78 ± 0.10. 



(4) 



This is shown in Figure [T] as the solid line; within errors it is essentially identical to equation ([2]) 
for —25 < My < —23, the interval over which we will be extrapolating the M, — L relationship to 
the most luminous galaxies in the sample. 

Figure [2] shows M,(L) based on a combination of the three relationships presented in Figured] 
plotted against M,(a") from equation ([T|) for the sample. The error bars along the M,(L) axis reflect 
the minimum and maximum predictions of M, given by the three relationships shown in Figure [H 
the central values plotted are the mean of the minimum and maximum predicted M,. The L and 
a predictors diverge at large L, with all three M, — L relationships predicting M, ~ 10^*^ Mq for 
the most luminous galaxies, while equation ([T]) predicts no values of M, larger than ~ 3 X 10^ Mq. 
The errors in M,(L) increase somewhat with galaxy luminosity but are much smaller than the 
differences between M,(L) and M,(cr), which approach an order of magnitude for some of the most 
luminous galaxies 



The differences between M,{L) and M,((t) cannot be reconciled by the IWyithd (120061 ) log- 
quadratic M, — a relationship. The asymmetric error bars in t he a-bas e d pre dictions of M, shown 
in Figure [2] r eflect the implied chan ge in predicted M, if the IWyithd (120061^ relation ship is used 
instead of the lTremaine et al.l (120021 ) log-linear M, — a relationship. The lWyithd (120061 ) relationship 
predicts slightly larger M, only for the largest a values (~ 30%), but still does not match the even 



^The error bars in Figure[2]do not include the systematic errors associated with the uncertainties in the individual 
relationships themselves. 
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larger M,{L) for the same galaxies. As expected, M,{L) and M,((j) do agree on average for the 
sample galaxies that actually have direct M, determinations, since it was this subset of galaxies 
that defined the relationships in the first place. 

The disagreement of the two M, predictors for the larger set of galaxies lacking direct M, 
determinations can be traced to changes in the form of the L — a relationship as a function of 
galaxy luminosity. Figure [3] shows this relationship for the sample galaxies. The typical a value 
appears to level off for large L; indeed, there appears to be little relationship between a and L for 
galaxies with My < —22. While most of the galaxies in this luminosity range are BCGs, other bright 
ellipticals show the same behavior. Put simply, the high luminosities of BCGs and other similarly 
bright ellipticals are not matched by similarly large velocity dispersions. The M, — a relationship 
thus predicts unexceptional black hole masses for these exceptionally luminous galaxies. 

This "saturation" iri a at BCG luminosities was noted in the BCG velo city dispersion study o f 
Oegerle Hoessell (jl99ll ). but it appears only weakly in the SDSS study of lBernardi et al.l (|2003l ). 
We suggest that this may be due to the use of different BCG luminosities, based on the analysis of 
the SDSS magnitudes of BCGs presented in Appendix [Al For the core galaxies, we find L ~ a"^ , 
a much steeper re lationship than the classic L ^ a"^. Specifically, a symmetrical least-squares fit 



(jPress et al.lll992l ) to the 99 core galaxies with My < —21 and having a a value produces: 



My = -2.5 (6.5 ± 1.3) log(cj/250 km s~^) - 22.45 ± 0.18. 



(5) 



However, since the L — a relationship appears to be nonlinear, even this fit may not be the best 
approximation for the most luminous galaxies. This result also contrasts with the relationship 
measured for power-law galaxies alone. 



My 



-2.5 (2.6 ± 0.3) log(cj/150 km s"^) - 20.30 ± 0.10. 



(6) 



The distribution of points with M, measurements shows what appears to be a bias in the 
BH sample: galaxies with My ~ —22.5 with measured M, have a higher-than-average a than 
typi cal galaxies at this lum inosity — or conversely have low luminosities for their a values (see 
also iBernardi et al.l l2006d ). The 7 galaxies with measured M, at My ~ —22.5 have average 
a = 311 lb 25 km s~^, while equation ([5]) predicts only ~ 25 km s~^ at My ~ —2 2.5 in agreement 
with the average a at this luminosity for the SDSS sample (IBernardi et al.l 120031 ). If a is the best 
predictor of M,, then the black holes in these galaxies should be on average (314/250)'^ ~ 2.4x 
more massive than is typical for galaxies with My ~ —22.5. The M, — L relationship in turn would 
be biased at the high luminosity end, and the large black hole masses predicted from L shown in 
Figure[2]will be over-estimates. Conversely, if L is the better predictor of M,, then then the M, — a 
relationship would be biased to predict lower M, than would be correct. 

The possibility that the galaxies with measured M, are a biased sampling of the L — o" re- 
lationship is echoed in Figure El For M, > 10^ Mq, M,{L) is on average greater that M,{a) for 
galaxies in the present sample. Lowering M,{L) by the bias factor inferred above, or increasing 



- 10 - 



M,(cj) by a similar factor would bring the average predictions into excellent agreement, however. 
Note the galaxies with measured M, in Figure [21 are presently in excellent agreement, since these 
are the very systems used to define the M, — a and M, — L relationships. 

Figure [2] also shows, however, that the large M,(L) predicted for the most luminous galaxies 
still deviate from M,{a) by a much larger factor than this putative bias. The strong curvature 
in L — a relationship leads to the upward curvature in M,{L) versus M,{a) well in excess of the 
selection biases implied by FigureO Any luminosity-based predictor of M, calibrated for My > — 22 
would still predict M, in excess of the M, — a relationship for My < —22, since a for the brightest 
galaxies does not increase with luminosity. 



4. Core Structure as an Independent Witness of M, 
4.1. The Cusp Radius 



Resolving whether L or cr is the best predictor for M, for galaxies with My < —23 will only be 
possible when real M, determinations can be made in this luminosity regime. Lacking this, we can 
attempt to obtain preliminary information by considering whether the central structure of galaxies 
may provide an independent witness to M,. We characterize the physical scale of the core by the 
"cusp radius," r^, which is the radius at which the negative logarithmic-slope of a galaxy's surface 
brightness pr o file re aches a pre-specified value 7'. This measure of core size was first proposed by 
Car olio et al.l (|l997l ): we will discuss it in detail in Appendix [Cl The core is also characterized by 
the cusp brightness, /-y, the local surface brightness at (//^ is expressed in magnitude units). 
In terms of the Nuker-law parameters, for 7 < 7' < /?, 



n 



7 -7 
/3-7' 



1/0 



(7) 



is then found directly from the fitted Nuker-law, 



= 2('3-7)/"4 



rb 



1 + 



(7-/3)/« 



(8) 



Car olio et al.l (|l997l ) advocated use of with 7' = 1/2 as a core scale-parameter. We show 



in Appendix O that using with 7' = 1/2, indeed gives tighter correlations with other galaxy 
parameters than the choice of rh as a scale parameter. While the Nuker-law rh is still used to 
calculate r-^ , we no longer use it directly as a measure of core size, in contrast to the analysis 



presented in lFaber et al.l (|l997l ). Lastly, we emphasize that since is generally well interior to it 
is not meant to describe the actual complete extent of the core; it is just a convenient representative 
scale. 
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4.2. Core Structure and Galaxy Parameters 



It has long been k nown that the physical scale of cores in early-type galaxies is correlated 



with galaxy luminosity ( Lauer 1985 : Kormendyl [l985 ) . This rel a,tionship may be d ue to the action 



of central black holes on the central distribution of stars (e.g., iFaber et al.l 119971 ). Figures H] and 



[5] show the relationships between u, L, and cusp radius, r^, for the present sample. The — a 
relationship is particularly steep, as varies by over two orders of magnitude, while a changes by 
only a factor of two. For core galaxies with My < —21, a symmetrical least-squares fit gives 



log(r^/pc) = (8.4 ± 1.6) log(o-/250 km s"^) + 1.99 ± 0.09, 



while the — L relationship is 



log(r^/pc) = (1.32 ± 0.n){-Mv - 23)/2.5 + 2.28 ± 0.04. 



(9) 



(10) 



Of the two relationships, L is the better predictor of r^, with only 0.31 rms scatter in logr-y, while 
the scatter of logr-y versus a is 0.63. Note that BCGs and non-BCG core galaxies appear to follow 
the same relationships between and a or L. For the non-BCG core galaxies. 



log(r^/pc) = (7.4 ± 1.2) log(<T/250 km s~^) + 2.00 ± 0.07, 



and 



log(r^/pc) = (1.94 ± 0.29)(-Mv - 22)/2.5 + 1.77 ± 0.06, 



while for the BCGs, 



and 



log(r^/pc) = (15.2 ± 7.5) log(o-/300 km s"^) + 2.80 ± 0.31, 



log(r^/pc) = (1.24 ± 0.17)(-Mv - 23)/2.5 + 2.28 ± 0.05. 



(11) 
(12) 

(13) 
(14) 



While the slopes of the relationships are different from those for the entire sample of core galaxies, 
there is no significant difference between the relationships within the parameter ranges in which 
BCGs and non BCGs overlap. A critical result that is evident in Figure [5] is that while BCGs have 
larger cores than less luminous core galaxies, they are completely consistent with the larger total 
luminosity of BCGs. In contrast, there is essentially no correlation between and a for > 300 
pc, as is evident in Figure luminosity is a much better predictor of core size in BCGs than a. 

The core is characterized by a surface brightness as well as a physical scale, thus one could 
also explore the relationships between I-y and a or L, but as we show in Figure [6l and r-y are so 
closely related that they can be regarded as interchangeable. The fitted relationship between the 
two parameters for core galaxies with My < —21 is 



fi^ = 2.5 (1.05 ± 0.07) log(r^/100 pc) + 16.23 ± 0.10, 
where /U^ is in units of 1/— band magnitudes per square-arcsecond. 



(15) 
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Lastly, and can be combined to estimate the stellar mass of the core interior to the cusp 
as = 7rI^r'^{M/ Ly), again using the conversion between mass and light given in the context of 
equation Symmetrical fits give the relationships between and L ox a as; 

log(M^/M0) = (1.35 ± 0.13)(-Mv - 22)/2.5 + 9.17 ± 0.05, (16) 

and 

log(M^/M0) = (8.6 ± 1.0) log(cr/200km s"^) + 8.55 ± 0.13. (17) 



4.3. Core Scouring and Black Hole Mass 

4-3.1. Black Hole Mass and 

The existence of the — a and — L relationships implies empirical relationships between 
and M,, given the separate M, — a and M, — L relationships. By combining equation ([T]) with 
equation ([9|) we find ~ or more precisely, 

M.(cr) =^ log(r^/pc) = (2.1 ± 0.4) log(M./lO^M0) + 2.7 ± 0.2. (18) 

At the same time, we can also combine equations (jlOp and (jH to find ~ M}'^^^'^, or 

M,(L) =^ log(r^/pc) = (0.96 ± 0.09) log(M,/lO^M0) + 1.9 ± 0.1. (19) 

Equations (|18p and (|19|) are inconsistent. The conflict between M,(L) and M,(cr) leads in turn to 
contradictory predictions for how the physical scale of cores is related to black hole mass. 

Comparison of the observed — M, relationship to the two inferred relationships presented 
above may offer a path to determining which of M,{a) or M,{L) is more accurate for the most 
massive galaxies. In the "core-scouring" scenario, cores are created by the orbital decay of a 
massive binary black hole, which would be formed during the merging of two galaxies. As the 
merger progresses, black holes in the nuclei of the progenitor galaxies are brought to the center of 
the merged sy stem by dynamical frict i on. W hile the center of the merger may initially be highly 



concentrated (jMilosavlievic Sz MerrittI l200ll ). as is the case for power-law galaxies, central stars 
interacting with the binary black hole are ejected from the center as the binary hardens. The 
ejection of stars erodes the steep central stellar density profile, creating a shallow cusp, or break 
from the steeper profile that still p ersists at larger radii. A core is the region of the galaxy interior 



to the break (cf. iLauer et al.lll995l ). 



Under this hypothesis, the relationship between core scale and M, ought to be more funda- 
mental than either of the — a or — L relationships alone. The action of the black hole mass 
on stellar orbits at the galaxy center creates the core structure directly, and the r^ — a and r-y — L 
relationships are then merely consequences of the separate M, — a and M, — L relationships. Ac- 
cording to this logic, we would conclude that the larger cores of BCGs are evidence of higher BH 
masses. 
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A major caveat standing in the way of this conclusion is that core scouring may not lead directly 
to a clean relationship between and M,. The binary B H ejects a total mass of stars, Mp,-, that 
is expected to be p roportional to the total merged M, (jQuinlanl Il996l : iMilosavljevic &: Merritt 
200 ll : iMerrittI 120061 ) . However, the resu l tant would depend on the radii over which stars are 
ejected from the center. Further, iMerrittI (|2006l ) presents simulations that show that core formation 
should be a cumulative process. Cores formed in one merger event will be depleted even further 
in subsequent mergers, presumably leading to even larger increases in that would reflect not 
only the total BH mass, but the integrated merger history as well. Under this hypothesis, cores 
resulting from successive dry mergers would be abnormally large compared to their BH masses, 
potentially explaining the extra-large cores of BCGs, which are thought to be formed by such 
multiple mergers. Under scrutiny, however, this explanation seems difficult to support, since the 
cores of BCGs show no excess compared to the luminosities of their host galaxies, and it is this 
latter quantity that is probably the best indicator of the amount of dry merging that any massive 
elliptical has experienced. In other words, the core masses of BCGs galaxies are the same fixed 
fraction of their total light as in other galaxies, not some amplified value driven by multiple mergers. 
Thus, we seem to be driven back to the basic explanation that the larger cores of BCGs are due 
simply to larger BH masses. 

Can we use actual core data to identify the correct r-y — M, relationship? Figure [7] tries this 
by plotting versus M, for the 11 core galaxies for which there are direct determinations of M,. 
A symmetric fit to and M, for these galaxies has the form 



log(r^/pc) = (0.83 ± 0.25) log(M,/lO^M0) + 2.20 ± 0.10. 



(20) 



This equation is essentially consistent with equation (jl9p . the relationship inferred from M,(L), 
rather than equation (jlSp . which inferred from M,((t). At the same time, the scatter in Figure [7] is 
large, thus this result is sensitive to how the r-y — M, relationship is fitted. For example, if r^ is 
treated as the independent variable in an attempt to predict M,, given r^, then 



log(r^/pc) = (1.5 ± 0.8) log(M./lO^M0) + 2.20 ± 0.11 



(21) 



(although we express r^ as the dependent variable for comparison with the relationships above). The 
slope of this relationship is intermediate between that in equations (fTSl) and (fT9]) . For completeness, 
if M, is treated as the independent variable, which corresponds to the scouring hypothesis that M, 



determines r^, then 



log(r.y/pc) = (0.59 ± 0.18) log(M,/10^ Mq) + 2.19 ± 0.10 



(22) 



These three fits do not in fact suffice to identify the "correct" r-y — M, relation for four reasons: 
1) the various slopes differ considerably because the native scatter in the data is large; 2) w e are 
seeking the "true" underlying relationship (i.e., the "theorist's" question of lNovak et al.ll2006l ). but 
without a knowledge of cosmic scatter and its separate contribution to both M, and r-y, we cannot 
fit the data properly to find it; 3) the slopes in equations (|18p and (jl9[) were likewise meant to 
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embody "true" relations, but they were derived from prior fits tliat themselves suffered a similar 
ambiguity; and 4) the sample of core galaxies with measured M, is potentially biased in some way 
that is not understood (cf. Figures [3] to [6]) , and any new fit based on these galaxies might therefore 
not be representative. On this last point, we emphasize caution. While the galaxies with measured 
M, may on average have offsets in the parameter plots shown, this does not mean a priori that the 
directly observed — M, relationship is biased. The small number of core galaxies with measured 
M, plus the number of parameters in play means that understanding any biases must await a richer 
sample. 

Likewise, the sample of core galaxies with measured M, will have to be increased considerably 
before it can be used to convincingly discriminate between the M,(L) and M,{a) relations. Nev- 
ertheless, we may be able to obtain some guidance in advance of such observations by comparing 
M, estimated from to values estimated from L or a. Figure [8] shows the results of using either 
equation ()20p or (12ip to predict M, from r^, in analogy to Figure O which compared predictions of 
Af, based on a versus L. Both versions of the — M, relationship predict larger M, than would be 
inferred from M,((t). The symmetrically- fitted — M, in equation (I20|) appears to be consistent 
with M,(L), also predicting M, ~ 10^^ Mq for the most massive galaxies. 

Presently, the large scatter in the observed — M, relationship and the attendant uncertainties 
in any empirical relationship derived from it does not decisively favor M,(L) over M,{a). Equations 
(j20p and (j2ip however, on average predict greater M, than would be inferred from M,{a), while 
equation ()20p is consistent with the larger black hole masses implied by L for the most massive 
galaxies. At this early stage the — M, relationship thus may favor consistency with the M,(L) 
relationship. The fact that the scatter of on L is smaller than that on a (as would be expected 
if is produced directly by black hole scouring and M, correlates more closely with L) as well as 
the fact that M,(L) plausibly explains the large core of BCGs as being due to more massive black 
holes, whereas M,{a) seems to provide no ready explanation this, may offer additional support 
that M,{L) is more appropriate for the most massive galaxies. 



4.3.2. Black Hole Mass and Core Mass 

An alternative approach to explore the relationship between core structure and black hole 
mass is to compare the core mass, = irl-yr'^^M/ Ly), rather than r^, to M,. Although, as we 
noted earlier, and are closely related, so the relationship between M-y and M, will contain 
information similar to the — M, relationship, core mass should be a more direct indicator of the 
amount of core scouring and its relationship to black hole mass. If cores are created from power- 
law galaxies by core scouring following a dry merger, one might expect that the core mass would 
be approximate ly proportiona l to the black-hole mass. This conjecture is supported by N-body 



calculations by iMerrittI (j2006l ) , who argues that the core mass produced by scouring in a single 
merger is ~ /M,, where M, is the mass of the merged black hole and / ~ 0.5, largely independent 
of the mass ratio of the merging black holes; he also argues that the total core mass after dry 
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mergers should be given by / ~ 0.5A^. Direct estimation of the mass ejected from the core by 
scouring is much more difficult observationally than theoretically, because we do not know the 
state of the galaxy before the merger. Thus we will simply use as an "indicative" core mass, 
recognizing that the factor / relating indicative core mass to black-hole mass is very uncertain, but 
should be approximately independent of galaxy luminosity for core galaxies. 

Figure [9] shows the relationships between and M, as derived from the M, — a and M, — L 
relationships. By combining M,{a) (equation [1|) with the My — a relationship (equation [T7|) . we 
find 

M,{a) =^ log(M^/M0) = (2.2 ± 0.3) log(M./lO^M0) + 10.29 ± 0.18, (23) 
while the combinations of M,{L) (equation U]) with the — L relationship (equation I16p gives 

M,(L) =^ log(M^/M0) = (1.0 ± 0.1) log(M./10^MQ) + 9.39 ± 0.11. (24) 

The relation between indicative core mass and black-hole mass predicted by the M, — L relation 
is essentially linear, as expected, while the relation predicted by the M, — a relation is twice as 
steep. Moreover the ratio of indicative core mass to black-hole mass is ~ 2.4 at M, = IO^Mq 
from the M, — L relation, not far from the value of order unity that we might expect, while the 
corresponding value from the M, — a relation is ~ 19. It is difficult to devise dynamical models in 
which core scouring could be efficient enough to create cores with mass so much bigger than the 
black-hole mass. 



5. The Growth of the Most Massive Galaxies and Why M,{L) Might be Favored 

Over M,(cr) 

Having found suggestive but not conclusive arguments to prefer one relation over the other, we 
turn now to physical arguments for additional guidance. We stress again that the tension between 
M,(L) and M,{a) arises due to the breakdown, or curvature, in the L — a relationship at high 
galaxy masses. It is appropriate to inquire at least briefly into physical reasons why this might 
happen, and whether this offers insight into which of M,(L) and M,((t) might be preferred for the 
most massive galaxies. 

Curvature in a gravitational scaling relation may signal a breakdown in perfect homology in 
galaxy formation, which in turn may reflect a change in the relative importance of different physical 
processes as a function of galaxy size. One such effect, suggested some time ago, is an increase 
i n the importance of dissipatioriless ( i -e., "dry") mergin g in forming the most massive ellipticals 
(jBender et al.lll992l : iFaber et al.l 119971 : iNaab et al.ll2006l ). A second effect, following logically from 
the first, is a change in the nature of dissipationless mergers with galaxy mass. As hierarchical 
clustering proceeds, clusters of galaxies become more massive, and previously formed elliptical 
galaxies at the centers of these clusters merge. Each new round of merging thus increases the 
galaxy stellar mass along with the cluster dark halo mass in which it is embedded. The largest 
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ellipticals are thus produced by dry mergers at the centers of the largest clusters. These are BCGs. 
The above scenario is supported by the steep environmenta l dependence amo ng bright ellipticals, 



the brightest ones being found in the densest environments (jHogg et alj 12004 ) 



There are at least two trends that might contribute to a breakdown in perfect homology for dry 
merging to produce the observed curvature in the L — a relationship. The f irst is a change in the 



typica l orbital eccentricity of central merging pairs as cluster mass grows. iBovlan-Kolchin et al 



(j2006l ) have suggested that head-on collisions may become more frequent when massive clusters 
merge, and their N-body simulations indicate less loss of energy from stars to dark matter in such 
collisions. The resultant stellar remnants are puffed up in radius and have significantly lower stellar 
velocity dispersions compared to encounters with normal orbital geometry. A second effect, not 
considered by them, is the fact that the ratio of cluster velocity dispersion to internal galaxy velocity 
dispersion also rises as clustering proceeds. This app ears to happen because gas cooling is reduced 



in large dark-matter halos (e.g., iBirnboim &: Dekell 12003). which means that the baryonic masses 



of central galaxies grow more slowly than their dark-matter halos. This is why the stellar velocity 
dispersions within BCG galaxies in large clusters are so much lower than those of their surrounding 



clust ers, whereas the same is not true of ellipticals in small groups (cf. Figure 2 of lBlumenthal et al 



19841 ). The net result is that central merging pairs will approach each other at relatively higher 
speeds, with the potential of injecting more orbital kinetic energy into the final stellar remnant. 
This would also cause the remnant to puff up and have smaller final velocity dispersion. 

The relative importance of these two effects can only be decided using realistic two-component 
N-body simulations containing both stars and dark matter that are appropriately embedded in a 
cosmological clustering scenario. It seems probable that both effects will be found to play a role. 
The point for now is that there are at least two reasons to expect non-homology in dissipationless 
mergers, and thus two reasons for curvature in the L — a relationship. 

If this logic is correct, it points towards M,(L) as being the proper scaling law for massive 
galaxies. That is because the major growth in black hole mass during dissipationless merging occurs 
by merging black holes as the galaxies themselves merge. With little mass accretion directly onto 
black holes during this stage and no attendant star formation, black hole mass should increase in 
proportion to galaxy mass. The ratio of stellar mass to black hole mass is constant over dry merging, 
consistent with M, scaling linearly with galaxy mass. Conversely, for M,((t) to be maintained over 
dry merging, given the plateau in a at high galaxy luminosity, would require one of the two merging 
black holes to be ejected from the galaxy as a common occurrence. These arguments provide yet a 
another reason to prefer the M,{L) relation over M,((t). 

Regardless of which mechanism is dominant for determining the velocity dispersion of the 
merged galaxy, "puffing-up" of the remnant does appear to happen in the most luminous galaxies, 
supporting an overall scenario in which velocity dispersion in the largest galaxies does not increase 
stro ngly via mergers. The evidence for this comes from the effective radii of the largest galaxies (see 



also iBernardi et al.ll2006bl ). Figure [TO] shows the Re — L relationship for the whole sample, where Re 
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is the effective radius measured from r^/^-law fits to tho se galaxies that have ground-based surface 
photometry extending to large radii (jLauer et al.ll2007al ). For low galaxy luminosity, the the mean 
Re — L relationship is relatively shallow. We find 



log(i?e/pc) = (-0.50 ± Om){Mv + 21)/2.5 + 3.62 ± 0.04 



(25) 



based on fits to just power-law galaxies. This stands in contrast to the steeper relationship defined 
by core galaxies with My < — 21, 



log(i?e/pc) 



-1.18 ± 0.06)(Mv + 23)/2.5 + 4.27 ± 0.02 



(26) 



The transition between the two forms occurs at My ~ —22, which corr esponds to the luni inosity 
at which the average central structure changes from power-law to core (jLauer et al.ll2007al ). This 
is also the scale at which a starts to plateau in Figure [3] — the leveling-off of the L — a relationship 
is thus associated with a rapid increase in Rg with L not seen in less luminous galaxies. This is as 
predicted if extra energy is injected into these galaxies by merging: the objects will have lower a 
but larger i?e. 

It should be stressed that the arguments presented in support of the M,[L) relation in this 
section apply only to bright ellipticals, which are those produced by dissipationless merging. The 
mass-accretion processes that built black holes when galaxies were younger were drastically different 
and might have obeyed different scaling laws. The M,{a) law might be a better fit to such galaxies, 
which in general will be smaller than the objects considered in this paper. The broader point of this 
discussion, however, is that "non-homology" processes may have affected the growth of galaxies 
generally, with the result that a single black hole scaling law with global galaxy properties might 
not fit all galaxies. 



6. The Space Density of the Most Massive Black Holes 

The preceding sections have presented suggestive if not conclusive reasons to suspect that the 
M, — L relation might be a better predictor of black hole mass than the M, — a relation for the 
most massive galaxies. 

1 . The velocity dispersions of the most massive elliptical galaxies rises slowly if at all with galaxy 
luminosity implying that their black holes are no larger than those of much smaller ellipticals 
if M, — fj is the governing relation — this seems rather surprising. 

2. The core-scouring model says that should correlate directly with M, whereas correlations 
between and L or a should be secondary; if so, the smaller scatter of on L validates L 
as the more accurate predictor of M,. 

3. The M, — L relation offers a simple explanation for the large cores of BCGs in terms of bigger 
black holes whereas the M, — a relation offers no such ready explanation. 
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4. Either of the two fits of M, on for core galaxies with measured BHs (equations 1201 and I2ip 
predicts large M, when extrapolated to luminous galaxies, in better agreement with M, — L 
than with M, — a. 

5. The largest elliptical galaxies are believed to be formed by dry merging, which predicts that 
black hole mass should grow in proportion to stellar mass; the observed M, oc L relation is 
thus the simplest relation predicted on these grounds. By contrast, the saturation of black 
holes mass in the largest galaxies that is predicted by the M, — a relation requires that one 
of the two merging BHs be ejected from the galaxy as a common occurrence, which may not 
by natural. 

6. The largest elliptical galaxies are likely formed by dry mergers at the centers of massive 
clusters; non-homology merger arguments plausibly explain the low velocity dispersions (and 
large radii) of such galaxies, but the same arguments then imply that a ought not to be a 
fundamental parameter for predicting black hole mass in the biggest galaxies. 

7. As noted in the introduction, there is evidence from AGN physics that M, ~ lO^^M© in some 
systems. 

Although these lines of argument are not conclusive, they motivate us to consider the impli- 
cations for the local BH mass function should the M, — L relation prove correct. The differences 
from the M, — a relation are large, as we shall show. In this section we first compute these two 
mass functions based on M, — a and M, — L, and then we compare the resulting mass functions 
to estimates of the relic black hole mass function based on the space density of QSOs as a function 
of luminosity and redshift. 



6.1. The Black Hole Cumulative Mass Distribution Functions 



We first compute the black hole mass function by combining the M,(a) predictor with the 
velocity-dispersion function (the space density of galaxies as a function of velocity dispersion). We 
then repeat the calculation, but thei i using M,(L) combine d with the galaxy luminosity function. 
Our analysis follows the precepts of lYu &: Tremaind (|2002l ). departing in the choice of dispersion 
functions. Both calculations use the same formalism, thus for the sake of generality we denote the 
log of either the velocity dispersion, cr, or luminosity of the galaxy, L, by s and assume that the 
correlations of BH mass M, with either parameter can be formalized through the statement that 
the probability of finding a galaxy a given BH mass is 

[log(M.)-F(s)]2- 



dP{M,) = (27rA' 



-1/2 



exp 



2A2 



dlog(M.), 



(27) 



where F{s) is the ridge line of either log(M, 
is then 

dn{M,) ,^_^2n-i/2 



s relation. The number of BHs near a given mass 



dlog(M.) 



(2^A' 



dn 

— exp 
ds 



[log(M.) - Fjs)]' 
2A2 



ds, 



(28) 
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and the cumulative distribution is 

n(M,) 



A/, 



dn{M,) 
. dlog(M. 



-d\og{M, 



(29) 



For the dispersi on-based predictor, we s tart with the lSheth et alj (|2003l ) SDSS-based velocity- 
dispersion function. iBernardi et alJ (j2006al ) reproces sed the SPSS data and recovered a number 
of galaxies with larger dispersions than those used in ISheth et al.1 (|2003l ) . We use that set of high 
dispersion galaxies to compute a cumulative dispersion function above a = 350km s~^ by 



n(> a) = N{> a)/V, 



(30) 



where A^(> a) is the number of galaxies w ith dispersions greater than a and V is the Sloan 
survey volume given bv lBernardi et alj (|2006al ) as 3.34 x 10"'^Mpc (for Hq = 70 km s~^ Mpc"^ and 
z < 0.3). This cumulative function is well approximated by a power law. Differentiating it gives 
an estimate of the dispersion function above a = 350km s^^ of 



dn{M,) 

(ilog(cr) 



6.67 X 10" 



350km s 



-10.27 



Mpc" 



(31) 



Equation [3T] predicts about 10 times as many ga l axies with dispersions greater than 400km s~^ 
as the Schechter f unction fit given in ISheth et al.l ( 20031 ). Above a = 350km we add it to the 
Sheth et al.l (|2003l ) dispersion function in the analysis of equation (j28p . 



We think this is the best available estimate of dn/dlog{a) for earl y-type galaxies a t zero 
redshift. We combine the Sheth/Bernardi dispersion function with the iTremaine et al.l (|2002l ) 
M,(cj) predictor (equation [1]) F[s) = 8.19 -|- 4.02 x x where x = s — log(200 km s~ ^), in eq u ation s 
[28] and [29l As an alternative, we consider the a - BH mass function using the IWvithd (|2006l ) 
predictor F{s) = 8.11 -|- 4.2 x x + 1.6 x x^. We illustrate both cumulative BH mass functions 
so derived in Figure [TTl Choosing the Wyithe predictor only increases the predicted BH number 
density modestly near M, = 3 X lO^M©. 



I nclusion of the cosm ic scatter in the M, relationships is crucial (|Yu &: Tremainelbood : Iyu k Lul 



20041 : iTundo et al.ll2006l ). The total population of galaxies at any given L or a will host black holes 
with a range in M,. The final BH mass functions thus are not a simple "relabeling" of the L or a 
distributions with black hole mass, but are rather a convolution of the of these distributions with 
an assumed distribution of M, at constant L or a. This convolution is especially critical at a large 
BH mass, where both the galaxy luminosity and dispersion functions decline rapidly. As a result of 
cosmic scatter, most of the high mass BHs actually come from "modest" galaxies with unusually 
large BHs for their luminosities or dispersion s, as compared to t he expected contribution of massive 
black holes from the most massive galaxies (jLauer et al.ll2007bl ). 



We emphasize that A above is the scatter a bout the mean r elatio n due to cosmic scatter in 
the relation and not due to measurement error. ITremaine et al.l (|2002l ) notes that the A = 0.30 
scatter about their derivation of M,{a) might be entirely due to measurement error, leaving no 
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room for cosmic scatter. In plots of BH mass functions in Figure [TT] we illustrate illustrate 3 values 
of the cosmic scatter for the Tremaine and Wyithe results: A = 0, 0.15 and 0.30. This scatter is 
probably not larger than 0.30, but may be considerably smaller. 

We compute t he th e L-based black hole mass function by the same approach, starting with 
the lBlanton et al.l (120031 ) SDSS g' galaxy luminosity function. We convert th e q' galaxy luniinosit y 
function to band using g' = V ^ 0.41 — suitable for E galaxies at z = (Fukugita et al. lOod). 
Comparison of the Blanton et al. galaxy luminosity function with the iPostman &: Lauer 
BCG survey suggests that the Blanton work undercounts BCGs. We argue in Appendix Rl that 
this may be due to the effects of excessive sky subtraction on the most luminous galaxies. To 
correct for this, we have added an estimate of the s pace density of BCGs as a function of V-band 
luminosity to the Blanton et al. sample. We used the lPostman Laued (jl995l ) BCG sample, which 
is volume- limited, to construct an estimate of the luminosity function in Rc-, transforming it to V 
using Kc = y — 0.55. Using the combined dn/dMy, we determine the number of BHs above a 
specified mass M, f rom equations [291 and [28l where we set s = My and F(s) is the right-hand side of 
equation ^ — the lHaring Rixl (j2004l ) predictor. Because iBlanton et al.l (120031 ) represented their 
luminosity function as that observed at z = 0.1, rather than the present epoch, both k-corrections 
and co rrections for e v olutio nary fading are required. These two terms fortuitously cancel each 
other: iBlanton et al.l (|2003l ) show that their sample dims by 0.2 mag in g' from z = 0.1 to the 
present, while the filter k-correction to transform to z = is —0.20 mag. Lastly, we use A = 0.25 
and 0.50, larger values than were used for M,{a), given the larger scatter in M,{L). 

The cumulative BH mass functions based on the two different methods are shown in Figure 
im Near BH masses of 10^ Mq (luminosities near L*) and lower, the L-based function overpredicts 
BH numbers by a factor of two and larger. In part, this is due to the fact that the L-based mass 
function includes galaxies that are disk-dominated, thereby overestimating the numbers and masses 
of BH at lower masses, since the L — M, relationship is based on bulges and elliptical galaxies. At 
higher masses this correction is negligible since most galaxies are ellipticals or bulge dominated SOs. 
There is also the possibility that the L — M, relationship is biased in the sense that the galaxies 
with black hole mass determinations have larger velocity dispersions than the average values for 
galaxies of their luminosities (see Figured]). If so, this would cause the L-based mass function in 
Figure [TT] to shift to the left, in better agreement with the cr-based mass function. 

For M, > IO^Mq the L and a mass functions diverge. The L based mass function predicts 
a local density of the most massive black holes that is about an order of magnitude greater than 
would be inferred from the M, — a relationship for M, > 2 x IO^Mq. This disagreement was 
foretold in Figure [2] — the present Figure [TT] simply recasts the disagreement between M, — L and 
M, — 0" in terms of the BH mass function. 
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6.2. The Black Hole Distribution Function Inferred Prom QSOs 

There are two different approaches that can be used to infer the present BH mass function from 
quasar counts (specifically from t he joint distri bution of quasar numbers as a function of redshift and 



luminosity). One line, started bv lSoltanI (jl982l ). relies on energy conservation. Under that argument, 
the energy emitted by the quasars at any redshift propagates through the universe declining in co- 
moving density due to the redshift of the photons as (1 + z)~^, and thus behaves exactly like a 
background. Hence the observed quasar flux translates directly to the total energy emitted given a 
known redshift distribution of emitters, and further translates into the t otal mass accreted by blac k 



holes, given their radiative efficiency. An alternate approach, started by Small Sz BlandfordI ( 1992 ) 



is to assume that all quasars go through a phase in which they accrete at the Eddington limit, 
followed by a period of slower or intermittent accretion according to a universal model dependent 
on BH mass and time. This assumption, together with a continuity argument (the number of BHs 
at a given mass changes only by accretion and merging) permits the recovery, not merely of the 
local BH density, but also of the local BH mass function. This second approach achieves a more 
detailed result than the Soltan argument, but at the expense of additional assumptions. 

A third more limited approach, which shares some logic with Small-Blanford, is to assume 
simply that the BHs of known mass accrete near the Eddington limit for some period and to ignore 
their fainter growth period — the so called "lightbulb model." In this model quasars are either 
on or off. Because the number of luminous quasars in the universe varies strongly with time, the 
model doesn't count BHs directly, rather it counts those that are accreting. In what follows, we 
evaluate the lightbulb model at z = 2.5 where the top end of the quasar LF is largest. We assume 
that no BHs above W^Mq are destroyed by mergers since that time, so the fall-off in the LF is due 
to a halt in accretion. At earlier times the BH mass function may be evolving. The full-width half 
maximum of the bright quasar LF is about 10^ yrs. So our lightbulb model implicitly assumes that 
every massive BH accretes for 10^/ yrs, where / is defined as the duty fraction, and then shuts off. 
This idealization ignores low mass BHs and low-level accretion. Another significant limitation is 
that the model provides no procedure to identify the mass below which it fails, although that failure 
is implicit within the assumptions: some of the lower luminosity quasars must be high mass BHs 
accreting at less than their Eddington limit. A third limitation is that the model is fundamentally 
inconsistent: the assignment of mass corresponding to a quasar luminosity gives the instantaneous 
mass of an accreting BH, while the present day mass function depends on its final mass. If the 
quasar is "on" for less than the Salpeter time (for BH mass to e-fold in Eddington- limited accretion), 
then the problem is small, but if it is on for much longer than the Salpeter time it is catastrophic: 
the present BH mass may be much larger than that assigned to the quasar. 

Nonetheless, the lightbulb model permits a comparison of the quasar LF with the number 
density of the most massive BHs in the local universe, with the duty fraction as a free parameter, 
under the assumption that the brightest objects i n the quasar luminosity function are accreting near 



the Eddington limit. This approach was used by lRichstone et al.l (|l998l ) to make a crude estimate 



of the duty fraction of BH accretion. We perform a similar analysis here. We start with the 
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Richards et alj (1200511 lum inosity function at z = 2.5, where the bright quasar density is greatest. 



The lRichards et al.l (|2005l ) fit reports number densities per magnitude at an AB absolute magnitude 
at rest-frame A =1500A. We integrate their fitting function from a given luminosity to infinity to 
get a cumulative number of objects brighter than that luminosity, we apply a bolometric correction 
of 5 (jMarconi et al.ll2004l ). and then we deduce a mass from the Eddington limit. This procedure 
identifies the number density of BHs greater than a given mass accreting at a given redshift. We 
compute this cumulative density at redshift z = 2.5 where the density of bright quasars is greatest. 
We divide th i s resu lt by the duty fraction f. We ado pt / = 0.03 based on the extensive work by 
Steidel et al] (|2002l l and kdelberger h SteideTI (|2005l l. This yields the hue labeled "hghtbulb" in 
Figure dH 

An improvement on both the lightbulb approach and the Small-Blandford approach is to use 
a physical model for the accretionary evolution of the BHs. One such model is th e merger-induced 



accret ion model that has been explored in detail by ISpringel et al.l (j2005l ) and iRobertson et al 



(1200a). They simulate the merger of disk+bulge+BH galaxies containing gas using the GADGET 
code, treating the BH growth by computing the Eddington-limited Bondi accretion rate at their 
smallest resolution element. They compute the luminosity of the accreting BH from the accretion 
rate under reasonable as sumptions about the radiative efficiency. Their simulations permit the 
development of a model (jHopkins et al.ll2006l ) that predicts the X-ray background and the zero- 
redshift BH mass function from the quasar LF. We believe the Hopkins model is a profound advance 
over simpler analyses. While it might turn out that their model does not correspond in detail to 
the quasar phenomenon, the approach may have broader utility. We summarize the salient points 
of their model below. 



The lHopkins et al.l (|2006l ) simulations exhibit very complex behavior of luminosity as a function 
of time for a given galaxy merger, but the time spent above a given luminosity turns out to 
be a universal profile over a wide range of galaxy or merger parameters, provided it is scaled 
appropriately with the peak luminosity and relic BH mass of the merger. For their simulations, 
the lifetime near a given bolometric luminosity L can be parametrized as 



dV^ilogL = t^exp(-L/L^) 



(32) 



where the timescale tg (a crude quasar lifetime) and luminosity scale Lg are functions of the peak 
luminosity Lp as follows. 

= 0.2Lp, (33) 



and 



1.37 I - 



-0.11 



Gyr. 



(34) 



The final or relic BH mass has a one-to-one correspondence with the peak luminosity given by 

-0.11 

X MEddiLp), (35) 



M, = 1.24 ( -4^] 
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where M^ddiLp) is the mass with an Eddington luminosity of Lp. An ensemble of objects with the 
same Lp should have a luminosity distribution proportional to the dt in equation | 



Hopkins et alj (|2006l ) use the model of quasar lifetimes described above together with a log- 
normal distribution of quasar birth rate per unit time to match the quasar luminosity function. We 
use their parameterization 



n 



cr* V2vr 



exp 



[\og{Lp/L,)f 
2ai 



(36) 



where n is the number of quasars born per unit comoving volume per unit time. iHopkins et al 
(|2006l ) find a good fit to the X-ray and optical quasar luminosity functions with 



LJz) 



L*(0) exp(A;r) 



(37) 



where r is the dimensionless lookback time r = Hq J dt and the other parameters are presumed 
constant. In what follows we use their best fit model with (logL,,, k, logn,,,, a^:) = (9.94, 5.61, 
-6.29, 0.91), with in solar units and n* in comoving Mpc~^Myr~^. 

We can compute the present day density of quasar relics by integrating the quasar birthrate 
over time at any specific mass or Lp. Therefore the cumulative density of BHs above a given mass 
M, is 

n(M,) 



r — oo 


/ n—dz 




Jo dz 



Following iHopkins et al.l (j2006l ) we set z. 



dlogLp. (38) 
3. We plot the result of equation (p8]l in Figure [TTl 



An important feature of the Hopkins model is that owing to the exponential distribution 
of time above a given luminosity (equation [32]) the quasar spends only a fraction of its lifetime 
accreting near the Eddington rate. For example, a W^Mq relic BH had a peak luminosity, Lp of 
3.03 X W^'^Lq and spent the time te = tgloge x Ei(5/e) = 15Myr above a factor of 1/e of Lp, 
where Ei{x) = exjp{—u)/u du is the usual exponential integral. The Hopkins model guarantees 
that the BH will accrete enough mass, but not too much, over its lifetime. 

Figure [11] permits us to compare the lightbulb and Hopkins models with the two relic BH 
mass functions. The BH mass functions diverge at about IO^Mq. The dispersion-based predictors 
predict considerably fewer BH at above W^Mq than the luminosity-based predictors. They are 
not con s istent with the lightbulb model; consistency with the Hopkins model is possible with the 
Wvithd ()2006l ) form of M,{a), but with the assumption of more cosmic scatter in the M, — a 
relationship than is probably realistic. The luminosity-based mass function is barely consistent 
with the lightbulb model, but probably overpredicts the AGN density compared to the Hopkins 
model. We thus cannot make a clear determination between the dispersion-based and luminosity- 
based BH mass pre dictors by comp a ring z ero redshift BH demographics to quasar demographics; 
however, the linear [Tremaine et al.[ ()2002l ) M,((t) relationship is disfavored in all of the present 
models to explain the QSO population. 
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An important caveat is that our calculations have neglected the effects of dry merging on the 
most massive galaxies after the epoch of QSOs. Merging might produce high mass black holes 
as a relatively recent phenomenon, thereby helping to reconcile the estimates made from M,(L) 
with the QSO population. Anoth er caveat for bo th results is the possibility of super-Eddington 
accretion among the biggest BHs (jBegelmanl |2006| ) . If common, super-Eddington accretion would 
make it very hard to make any estimates of the mass function of relic BHs from quasar LFs. 



7. Conclusions 

The M, — a relationship has come to be the "gold standard" for predicting black hole masses 
from galaxy properties due to its small scatter and its implications for illuminating the co-formation 
of galaxies and their nuclear black holes. At first sight, the M, — L relationship might be dismissed 
as a simple consequence of the Faber-Jackson relationship. With M, ~ and L ~ a^, one would 
expect something like M, ~ L. The larger scatter in the M, — L relationship further suggests that 
the M, — a relationship is really more fundamental. But as galaxy luminosity increases, a levels 
off and the basic Faber-Jackson relationship does not appear to hold. At BCG luminosities there 
are no direct measurements of M, and M,(o") versus M,{L) present contradictory extrapolations. 
The contradiction essentially begs the question, do these exceptionally luminous galaxies have 
exceptionally massive black holes? The M, — L relationship answers this in the affirmative, while 
for the M, — a relationship to be correct we must accept the puzzling result that the black holes 
in BCGs have relatively modest masses. But this question leads to a broader issue, namely, is is 
there a significant population of black holes with M, approaching 10^^ Mq in the local universe? 

The best way to answer these questions is to attempt to weigh the BHs in BCGs. With 
the advent of LASER-guided adaptive optics-fed spectrographs on 10m class telescopes, it is now 
possible to do this. This paper may therefore be premature. However, given the high attention 
to the M, — a relationship and its implications for galaxy formation, we believe that advancing 
the implications of M, — L relationship offers an important alternative view that should not be 
overlooked. Lacking hard measures of M, in the most massive galaxies, we have marshalled a 
number of less-direct arguments that this hypothesis may be favored for the most luminous galaxies. 

The first set of arguments is based on the hypothesis that cores in the most luminous galaxies 
are created in a "core scouring" process in which a binary BH created in the merger of two galaxies 
evacuates stars from the center of the newly-merged product. There presently is little observational 
support for the creation of binary black holes in mergers, but abundant theoretical work shows 
that realistic cores can be created by binary black holes, and the prevalence of nuclear black holes 
in galaxies overall strongly argues that such binaries must be created as a natural consequence of 
mergers. If so, the physical scale of cores, which we have parameterized as may be an independent 
witness of M,, and thus use the large cores in BCGs as an indicator of their black hole masses. 

Based on central structural parameters derived from HST observations, we find that the large 
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cores in BCGs are commensurate with their high luminosities, while cr is a poor predictor of 
for > 300 pc. The scatter in the — L relationship is much smaller than that in the — a 
relationship, again implying that L and core scale are more closely related. The observed — M, 
relationship for 11 core galaxies with directly determined black hole masses has large scatter, but 
appears to be more consistent with the M, — L rather than the M, — a relationship. Lastly, the 
core masses in BCGs are over an order of magnitude larger than the black hole masses estimated 
from M,(cr), but are only a few times larger than those estimated from M,(L); making such large 
cores with the smaller cj-based black hole masses would be a strong challenge for the core-scouring 
hypothesis of core formation. 

The second set of arguments comes from considering theoretical arguments concerning whether 
or not L rather than a should predict M, in BCGs. The favored origin of BCGs is that they are the 
remnants of dissipationless purely-stellar mergers of less-luminous elliptical galaxies, augmented by 
ongoing galactic cannibalism of elliptical galaxies in the rich environments at the center of galaxy 
clusters. The plateau in the L — a relationship plus but the steeper L — relationship at high 
galaxy luminosity presented here strongly favor this formation scenario. The luminosity of a BCG 
is the sum of the luminosities of its progenitors. Similarly, setting aside the possible ejection of 
nuclear black holes in the final stages of a merger, the final nuclear black hole mass should be the 
sum of the progenitor black holes. Stated more directly, the ratio M,/L should be largely invariant 
over dissipationless mergers, leading to M, ~ L at the high end of the galaxy LF. In contrast, a 
appears to be nearly constant over these mergers. In effect, even if a relationship between M, and 
a were set up in the initial stages of galaxy formation, it might not survive in a dissipationless 
merging hierarchy. 

A final argument comes from attempting to infer the z = space density of the remnant black 
holes associated with the most luminous QSOs seen at z ~ 2. As noted in the Introduction, the 
properties of the broad-line regions in the most luminous QSOs argue that they are powered by 
black holes with M, ~ IO^^Mq. The heating of the intra-cluster medium in the richest galaxy 
clusters may also demand that some black holes in BCGs approach this mass. The critical part of 
this analysis is understanding how to correct the QSO space density for QSO luminosity evolution. 
The remnant black holes last forever, but the QSOs represent only those BHs made visible during 
an epoch of high mass- accretion, which presu mably lasts only a small fraction of the age of the 



universe. We used the iHopkins et al.l (120061 ) simulations to estimate the QSO lifetimes. The 
resulting shape and space density of the bright end of the QSO LF falls between the higher space 
density of the most massive black holes implied by M,{L) and those implied by M,((t), while the 
simple "lightbulb" model of QSO duty cycles favors the M,(L) relation. This treatment is sensitive 
to the assumed amount of cosmic scatter in both M, relationships; however, it appears difficult for 
the log-linear M, — a relationship to explain the the observed space densities of the most luminous 
QSOs without assuming that its cosmic scatter is larger than is likely to be the case. 
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A. The Luminosities of BCGs and Comparison to SDSS Magnitudes 



Our analysis depends critically on the accuracy of the absolute luminosities of the brightest 
galaxies in the sample, such as BCGs. This is underscored by the bright-end disagreements of our 
L — a relation ship and galaxy lurn inosi ty function with tho se based on Sloan Digital Sky Survey 
(SDSS) data (IBernardi et alJ l2003l and iBlanton et alJ l2003l . respectively). We thus describe the 
derivation of our BCG luminosities, and compare them to magnitudes based on the SDSS for 
BCGs in common. We conclude that the SDSS BCG magnitudes are strongly affected by sky 
subtraction errors, causing them to be biassed to significantly fainter values. 



The BCGs in the present sample come from the 



Laine et al 



mmHST BCG study. This 



program, in turn, was based on the volume-limited iPostman Sz Laueij (|l995l ) BCG sample, which 
provides ground-based profiles and aperture photometryj^ As outlined in ^2.11 we derive apparent 
magnitudes of the BCGs by fitting the classic r^/'* form to the inner portions of the brightness 
profiles, where the inner limit of the fit was set to avoid seeing and the ou ter limit was specified 
to avoid portions of the profile that appeared to rise above the r^/^ fit. Graham et al.l ( 19961 ) 
showed that the BCG profiles could be described by single-component SersicI ( 1968 ) forms, but 
ones that often had index n > 4. The apparent magnitudes, which were derived by integrating the 
r^/^-law over radius, thus if anything are underestimates of the total BCG fiuxes. An alternative 



to thi s procedure would be to integrate the Sersic forms, however, as is shown in iGraham et al 



(|l996l ). the Sersic Re and n values are closely coupled, such that large n is matched with large Re, 
The implied total magnitude strongly diverges as n increases, and must essentially be regarded as 
unphysical extrapolations because the derived Re is typically well outside the actual radial domain 
of the surface brightness profile for large n; this is not true for n = 4. A contrasting treatment that 
occurs in much of the literature is based on the presumption that BCG must be completely well 
described by r^/'^-laws (in contrast to other giant elliptical galaxies, which also have n > 4), and 
that Sersic n > 4 is really the signature of an intracluster light component that must be subtracted. 
We conclude that an unambiguous procedure to derive total BCG luminosity does not presently 
exist. Our procedure of deriving magnitudes from just the inner portion of the profile that is well 
described by an r^/^-law again should give a sensible lower limit to the total luminosity. 



The high luminosity end of the 



Blanton et al 



(l2003l'l lu minosit y function falls well be low the 



BCG space densities measured by iPostman fc Laueii ()l995l ). The iBernardi et al.l (120031 ) L — a 
relationship shows no plateau at its bright limit. These discrepancies would both be explained if 
the total magnitudes for BCGs in the SDSS database are signifi cantly underestimated. We checked 
this hypothesis by examining the 25 IPostman &: Laueii (^19951 ) BCGs present in the SP SS DR4 
databases In detail, we compared the total Rc apparent magnitude in lLaine et al.l (|2002l ) against 



Postman fc Laued (|l995l ) did not actually pu blish their BC G surface brightness profiles, but they were presented 



graphically in the BCG profile analysis of Graham et al. ( 19961 ). 

®Two additional BCGs were mis-identified in the SDSS database as stars. 
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the SDSS r "model magnitude," (which in almost all cases is the most luminous total magnitude 
provided by the SDSS database) transformed by Rc = r — 0.31. The results are shown in Figure 
[T2] as a function of effective radius (based on our fits). A strong systematic trend is evident 
such that larger galaxies have greater offsets between the two total magnitudes. The median 
r - 0.31 - Rc value is 0.54 mag and rises to 1.57 mag for the NGC 6166 (the BCG in A2199). As 
an additional check, we also compared the SPSS r magnitudes against the maximum Rc aperture 
magnitude published by lPostman &: Laued ()1995l ). The maximum aperture radius was not defined 
in any rigorous way and does not correspond to any fixed fraction of the total galaxy flux, but 
the magnitude is a model- independent integration of all the flux within the published radius. The 
median difference between the SDSS model r magnitudes (transformed to Rc) and the maximum 
aperture magnitude is 0.24 mag and rises to values over a full magnitude for the largest galaxies. 
This demonstrates directly that the SDSS model magnitudes for the galaxies in question cannot be 
regarded as total magnitudes. 

Conversations with a number of experienced users of the SDSS database for bright galaxies 
warned us that the SDSS pipeline measure d sky levels on angular sc ales too small to accommodate 
bright nearby BCGs of the sort observed bv lPostman Sz Laueii (jl995l ). and indeed the results shown 
in Figure [12] strongly suggest that a sky-subtraction error affects the SDSS magnit udes. As a check, 
we plot the SDSS r surface brightness profiles against the lPostman &: Laueii (jl995l ) profiles for three 
of the BCGs with the largest ma gnitude differences i n Fig ure [T3l The SDSS profiles agree well at 
small radii but all fall below the iPostman &: Laueii (|l995l ) profiles at large radii, consistent with 
large sky subtraction errors. 

T he large SDSS sky-sub trac tion errors for bright galaxies may have important implications 
for the lBernardi et al.l (|2003l ) and lBlanton et al.l (|2003l ) studies, but exactly how is not clear. Both 
SDSS studies are based on galaxy samples with higher mean redshifts than the iPostman &: Lauer 
(I1995I ) sample. Their BCGs should be angularly smaller and thus be less vulnera. b le to sky sub- 
traction errors. Typical BCGs in the SDSS sample are listed in the iMiller et al.l (|2005l ) sample 
of galaxy clusters identified from SDSS galaxy catalogues. Figur e [T4l shows a histog ram of SDSS 
model r magnitud es (converted to My) for the BCGs identified by lMiller et al.l (120051 ) compared to 
a histogram of all IPostman &: Laued (|l995l ) BCGs with My based on their tota l Rc magnitudes. 
There is a clear offset between the two samples, with the iPostrnan &: Laued (119951 ) BCGs appearing 



to be typically one magnitude b righter than the Miller et al. " (j2005l ) BCGs. However, a histogram 



of SDSS r magnitudes for th e 25 Postman Sz Lauei] ( 19951 ) BCGs in common with SDSS agrees well 



with the iMiller et al.l (|2005l ) BCG histogram, yet these are the magnitudes shown to in be error. 
We conclude that the total magnitudes of nearby SDSS BCGs are wrong. 



B. The Use of Catalogued 2MASS XSC Apparent Luminosities of BCGs 



After the first version of this paper was posted on astro-ph. iBatcheldor et al. (2006) presented 
a M, — L relationship for BCGs based on apparent magnitudes extracted from the 2MASS Ex- 
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tended Source Catalogue (XSC) (jjarrett et al.ll200Cll . 120031 ) . The implied NIR luminosity differential 
between BCGs and other giant elliptical galaxies is greatly reduced from that of the present work. 
As a result, the plateau in the L — a relationship presente d here is greatly r e duced in the NIR and 
the conflict between M,{a) and M,{L) is thus resolved. iBatcheldor et al.l (j2006l ) further suggest 
that the relatively higher luminosities inferred from the optical photometry may imply that the 
envelopes of the BCGs are extremely blue. 

We have not conducted a complete comparison of the present photometry with that pro- 
vided by the 2MASS XSC, but a spot check of a few systems makes it clear that the 2MASS 
i magery from which the cat alogue magnitudes were derived is extremely shallow compared to that 
of iPostman &: Laueii (|l995l ). which is the source of the R band optical photometry (transformed 
to V) used in thi s pape r. The most likely explanation for the difference between the present and 
Batcheldor et al.l (j2006l ) results is that the 2MASS images are simply not deep enough to obtain 
accurate total luminosities of the BCGs, at least as represented by the automatic reductions used 
to generate the XSC magnitudes. 

Figure 15] shows the J band 2MASS image of NGC 2832, the BCG in A0779 compared to 
the central portion of the R band image obtained bv IPostman &: Laueii (|l995l ). The J sky level is 
15.67 magnitudes per square arcsecond versus the R band sky level of 20.90. Accounting for the 
R — J = 1.68 color in the center of the galaxy implies that the J band has a sky level effectively 26 x 
brighter. Further, the 2MASS image is a 7.8s exposu re obtained with a 1.3m telescope as compared 
to the 200s R image obtained with a 2.1m telescope ([Postman Laueiill995l ). The J image is thus 
considerably shallower than the R image as is evident in Figure [151 The galaxy envelope in the J 
band image disappears into the noise at radii at which it is still clearly present in R. This problem 
is exacerbated in the H and K bands, which have yet brighter sky levels. 

Despite the much brighter J sky level, the J and R band profiles of th e A0779 BCG are 
consis tent, as is evident in Figure [TBI which compares the R profiles measured by [Postman &: Lauer 
(|l995l ) to J profiles derived by us from the 2MASS archive images for the three BCGs shown earlier 
in Figure [T3l The final J band isophote shown for A0779 falls fully ~ 6.8 magnitudes below the sky 
level, but is still in agreement with the R band profile within the large error bars, which represent 
a 0.004 magnitude error in the 2MASS J sky level. 

An r^/^ law fitted to the J band profile of A0779 yields rrij = 9.04, only 0.08 magnitudes fainter 
than mj = 8.96, e stimated by subtra c ting t he R — J = 1.68 color of the central isophotes from niFt 
obtained from the [Postman &: Laueii (|l995l ) photometry. These values are in poor agreement with 
the XSC isophotal (</mfc2o/e)' extrapolated {Jniext) "^J values of 9.78 and 9.67, respectively, 
however. The isophotal radius, rk20fe, is 53 "3, well interior to the limits of the surface photometry 
shown in Figure [T6j The extrapolated magnitude is based on a Sersic fit to a surface brightness 
profile generated by the XS C pipeline. However , even for giant elliptical galaxies the XSC Sersic 



index is limited to n < 1.5 ([Jarrett et al.[ [20031 : Jarrett, private communication). The XSC cal- 



culation of total magnitudes thus assumes that the galaxies essentially have exponential profiles. 
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rather than the r^/^ form standard for giant ehi ptical galaxies. The XSC pipehne gives n = 1 . 17 fo r 



A0779, for example, while iGraham et alj (|l996l ) find n > 10 based on the lPostman &: Laueij (|l995l ) 
photometry. An exponential cutoff explains both the small difference between the XSC isophotal 
and total magnitudes, as well as the large deficit of both magnitudes compared to a total magnitude 
estimated from an r^l'^ law. 

A similar pattern is seen for the two other BCGs shown in Figure [161 For A2052, the BCG 
has Jrrn.2o/e — 10.92, corresponding to rfc2o/e = 38''4, and Jme^t — 10.60, while an r^/^ fit to the 
surface photometry recovered from the 2MASS archive image gives mj = 9.58, a full magnitude 
brighter and only 0.21 magnitudes dimmer than the mj = 9.36 inferred from the R band imagery 
with R - J = 1.79. For A2199, the BCG (NGC 6166) has Jm^ao/e = lO-^l, corresponding to 
i'k20fe = 50''0, and Jm.^^ = 1P;41, g enerated from a Sersic fit with n = 1.18 (Jarrett, private 
communication) : ICraham et al.l ( 19961 ) find n = 6.9. An r^/^ fit to the surface photometry recovered 
from the 2MASS archive image gives mj = 9.66, 0.75 magnitudes brighter, but 0.41 magnitudes 
dimmer than nij = 9.20 inferred from the R band imagery with R — J = 1.82. 

To summarize, our fits to the 2MASS J images for the three BCGs shown give mj values 
markedly brighte r than the XSC apparent magnitudes, but that are much more consistent with the 
R photometry of iPostman &: Laued (jl995l ). As a check, we found that the aperture photometry in 
Postman Sz Laueir ( 1995 ) interpolated to rk20fe was in excellent agreement with the Jmfcao/e values, 
assuming the R — J values derived by comparing the inner portions of the surface photometry 
profiles; however, rfc20/e is well interior to the limits of the optical photometry, thus the XSC 
isophotal magnitudes cannot be regarded as a total apparent magnitudes. The nearly exponential 
profile form used by the XSC pipeline explains why Jm^^^t is only modestly brighter than Jm^ao/e 
for the three galaxies shown, but is much dimmer than our fits, which assume Sersic n = 4, the 
standard value for giant elliptical galaxies. 

This analysis thus leads us to conclude that th e extremely blue BCG envelopes and the rel- 
atively modest NIR BCG luminosities advanced by iBatcheldor et al.l (|2006l ) are artifacts. We do 
note that the J surface photometry does fall below the R photometry at large radii in all three 
BCGs, but this occurs for J isophotes well within the noise and well below the even the errors in 
the sky levels. Any color gradients implied by the profiles presented in Figure [16] are not significant, 
and in any case are considerably smaller than would be implied by direct comparison of the XSC 
apparent magnitudes to the parameters given in this paper. 



C. The Selection of over r;, as the Core Scale 



The break-radius parameter in the Nuker-law has been used directly as the indicator of physical 
core scale in earlier work by our group. iFaber et al.l (|l997l ) showed that r;, is correlated with My 
for core galaxies. The scatter about the — My relationship is large, however. With the present 
much larger sample of galaxies, we were able to conduct extensive searches for other parameters 
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that might reduce the scatter, with the goal of better understanding how core structure varies with 
galaxy properties. 

Several plots and parameter fits were conducted to search for correlations between residuals 
of the rb — L and rj, — a relationships with the Nuker profile shape parameters a, (3, and 7. We 
also tried local values of 7, 7', measured at constant fraction of r^, interior to {^'{fff,) with 
/ < 1). A correlation emerged between the Nuker-law outer profile slope /? and the residuals in 
the r^ — L relationship, as is shown in Figure [T71 Larger than average cores correspond to galaxies 
with large (5 and vice versa. While one might be tempted to use this relationship to yield some sort 
of "/3-corrected" break radii, the form of the correlation suggests a simpler, more useful approach. 
In the Nuker law, marks the maximum curvature of the profile in logarithmic coordinates, but 
also the location where the profile slope reaches the average of 7 and (3. Since the range of 7 for 
core galaxies is greatly restricted compared to that of [3, this effectively means that rh will fall at 
steeper portions of the profile for steeper /?. This suggests that some sort of "cusp radius," r^, a 
radius at which the profile reaches a nominal slope value, 7', would provide a bett er description 



of the core scale than r^. The exact definition of we use is given in equation ([8]). ICarollo et al 



(jl997l ) had already suggested use of with 7' = 1/2, the choice that we adopt here. 



Figure [TSl compares the — L relationship for 7' = 1/2 (shown earlier in Figure [5|) and the 
more traditional — L relationship. A symmetrical fit to the latter yields 

log(rb) = (-0.71 ± 0.04) (My + 21.64) /2.5 + (2.05 ± 0.02), (CI) 



which can be compared to the — L relation given in equation (jlOp . For most galaxies /3 > 1, so 
in general, < r;,; on average ~ O.Sr^ for 7' = 1/2. We measured the rms scatter in the — L 
relationship for core galaxies with My < —21 as a function of 7' over the range 0.4 < 7' < 0.7. A 
broad minimum in the scatter of 0.31 in logr^ occurs at 7' ~ 1/2. This is significantly smaller than 
the 0.38 scatter in logr;, for the rb — L relationship. The reduced scatter in the — L relationship 
as compared to the rb — L relationship is clearly evident in Figure [I8l We have thus chosen to use 
over as the core scale. 

Evaluation of at 7' = 1/2 also leads to a clean separation of core and power-law galaxies. 
Since power-law galaxies have 7' > 0.5 at the HST resolution limit, they are naturally excluded 
from the analysis. The upper limits for for power-law galaxies however have the same physical 
values as their rb limits; since < rb for core galaxies, this may create a false impression that the 
upper limits of for power-law galaxies are more in line with the typical values of core galaxies. 

The issue of whether or not intermediate galaxies can be included within the class of core 
galaxies, or should be treated separately, unfortunately depends on which relationship is being 
considered. As is evident in Figure [TSl the intermediate galaxies with My < —21 appear to be 
evenly distributed about the the mean rb—L relationship; their mean residual about the relationship 
is — 0.06ib0.15 in log(rb). In contrast, the same galaxies fall preferentially to the compact side of the 
r^ — L relationship, now having a mean log(r^) residual of — 0.60ib0.11; given this systematic offset, 
we conclude that the intermediate galaxies should be treated separately from the core galaxies. 
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Table 1. Galaxy Parameters 
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NGC 0404 


SO- 


-17.19 


38 


\ 


0.23 


14.62 




NGC 0474 


SO 


-20.12 


164 


\ 


1.15 


15.20 




NGC 0507 


so 


-23.02 


311 


n 


2.22 


16.62 




NGC 0524 


S0+ 


-21.85 


253 


A 


1.57 


15.24 




NGC 0545 


BCG 


-22.98 


242 


n 


2.16 


16.52 


A0194-M1 


NGC 0584 


E 


-21.38 


207 


n 


0.95 


14.06 




NGC 0596 


E 


-20.90 


152 


\ 


0.63 


13.99 




NGC 0720 


E 


-22.20 


242 


n 


2.54 


17.22 




NGC 0741 


E 


-23.27 


291 


n 


2.46 


17.48 




NGC 0821 


E 


-21.71 


200 


A 


0.66 


13.75 




NGC 0910 


BCG 


-22.79 


257 


n 


2.21 


16.96 


A0347-M1 


NGC 1016 
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-22.90 


294 
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2.25 


17.01 




NGC 1023 


SO- 


-20.53 
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0.36 


12.90 




NGC 1052 
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-21.17 


208 
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1.46 


14.35 
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-20.13 
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14.09 




NGC 1316 
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-23.32 
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1.54 
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SO- 
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NGC 1374 


E 


-20.57 
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14.57 




NGC 1399 
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-22.07 
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16.76 
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SO- 


-21.37 
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1.47 


15.25 




NGC 1426 


E 


-20.78 
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0.71 


14.28 




NGC 1427 
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-20.79 
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0.61 


14.11 




NGC 1439 
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-20.82 
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0.71 


13.85 




NGC 1500 
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-22.75 
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1.99 


16.34 


A3193-M1 


NGC 1553 


SO 


-22.06 
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1.01 


13.54 




NGC 1600 
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-23.02 


337 
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2.82 


18.17 




NGC 1700 
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-21.95 
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1.01 
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NGC 2300 
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-21.74 
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2.12 


16.82 




NGC 2434 
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-21.33 
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0.64 


14.60 




NGC 2549 
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-19.17 
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13.98 




NGC 2592 
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-20.01 


265 
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0.82 


13.76 
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Table 1 — Continued 
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NGC 2634 
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-20.83 
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0.93 


14.57 




NGC 2636 
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-19.19 


69 
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0.87 


15.22 




NGC 2685 


S0+ 


-19.72 
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0.84 


14.16 




NGC 2699 
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-20.25 
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0.84 


14.14 




NGC 2778 
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-18.75 


162 
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0.67 


13.97 




NGC 2832 


BCG 


-23.76 


335 
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2.52 


17.11 


A0779-M1 


NGC 2841 


Sb 


-20.57 


206 
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1.09 


14.54 




NGC 2872 
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-21.62 


285 
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1.06 


13.65 




NGC 2902 


SO 


-20.59 
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2.15 
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NGC 2907 
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-21.23 
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1.22 


13.47 
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SO 


-19.73 
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NGC 2974 
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-21.09 
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13.77 




NGC 2986 
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-22.32 
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2.07 
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NGC 3056 
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-18.98 




\ 


0.80 


14.10 




NGC 3065 
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13.93 
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-21.95 
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NGC 3384 


so- 
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92 


\ 


0.65 
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NGC 3608 
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-21.12 
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1.31 
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NGC 3610 
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-20.96 
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0.64 


12.86 
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NGC 3613 
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-21. 


.59 


210 


n 


1. 


.65 


15. 


,65 




NGC 3640 


E 


-21. 


.96 


182 


n 


1. 


.47 


15. 


,39 




NGC 3706 


SO- 


-22, 


,31 


270 


n 


1, 


,60 


14. 


,16 




NGC 3842 


BCG 


-23, 


,18 


314 


n 


2, 


,48 


17. 


,40 


A1367-M1 


NGC 3900 


S0+ 


-20, 


,80 
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1, 


.16 


14. 


,25 




NGC 3945 


S0+ 


-20. 


.25 


174 


\ 


0. 


.59 


14. 


,19 




NGC 4026 


so 


-19. 


.79 
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\ 


0. 


.48 


12. 


,96 




NGC 4073 
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-23, 
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278 
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,13 


16. 


,55 




NGC 4121 
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.53 
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14. 
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NGC 4128 


SO 
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13. 
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NGC 4143 


so 
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\ 
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.88 


13. 
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NGC 4150 


so 


-18. 
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0. 


.85 


13. 


,87 




NGC 4168 
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-21. 


.80 


184 


n 


2. 


.26 


17. 
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NGC 4239 
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,50 
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1, 


,06 


16. 


,82 




NGC 4261 
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-22, 


.26 


309 
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2, 


.31 


16. 


,09 




NGC 4278 
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-21, 


,05 


238 
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1, 


,77 


15. 


,82 




NGC 4291 
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-20, 


.64 
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1, 


.63 


15. 
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NGC 4365 
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-22. 
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256 
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2. 


.15 


16. 


,53 




NGC 4374 
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-22. 


.28 
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n 


2. 


.11 


15. 


,67 




NGC 4382 


S0+ 


-21, 


,96 


179 


n 


1, 


,69 


15. 


,34 




NGC 4387 


E 


-19, 


.25 


104 


\ 


0, 


,54 


15. 


,13 




NGC 4406 


E 


-22, 


.46 


235 


n 


1, 


,90 


16. 


,03 




NGC 4417 


SO 


-18. 


.94 


131 


\ 


0. 


.94 


13. 


,96 




NGC 4434 


E 


-19. 


.19 


120 


\ 


0. 


.54 


14. 


,44 




NGC 4458 


E 


-19. 


.27 


103 


n 


0. 


.80 


13. 


,57 




NGC 4464 


E 


-18, 


.82 


127 


\ 


0, 


.54 


13. 


,92 




NGC 4467 


E 


-17, 


,51 


68 


\ 


0, 


.54 


15. 


,07 




NGC 4472 


E 


-22, 


,93 


291 


n 


2, 


,25 


16. 


,53 




NGC 4473 


E 


-21. 


.16 


179 


n 


1. 


.73 


15. 


,40 




NGC 4474 


SO 


-18. 


.42 


87 


\ 


0. 


.72 


14. 


,74 




NGC 4478 


E 


-19. 


.89 


138 


n 


1. 


.32 


15. 


,50 




NGC 4482 


E 


-18. 


.87 


26 


A 


2. 


.05 


19. 


,52 
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Table 1 — Continued 



Galaxy 


Morph 


(J 

Mv (km/s) 


P 


log(pc) 


(V-Band) 


Alt-ID 


NGC 4486 


E 


-22, 


.71 


332 


n 


2, 


.65 


17. 


,25 




NGC 4486B 


cE 


-17, 


.98 


170 


n 


1 


.08 


14. 


,44 




NGC 4494 


E 


-21 


.50 


150 


\ 





.54 


13. 


,40 




NGC 4503 


SO- 


-19 


.57 


111 


\ 





.63 


14. 


,42 




NGC 4551 


E 


-19, 


.37 


108 


\ 





.54 


14. 


,86 




NGC 4552 


E 


-21, 


.65 


253 


n 


1, 


.60 


15. 


,17 




NGC 4564 


E 


-20, 


.26 


157 


\ 


0, 


.63 


13. 


,43 




NGC 4589 


E 


-21 


.35 


224 


n 


1 


.40 


15. 


,41 




NGC 4621 


E 


-21 


.74 


225 


\ 





.54 


12. 


,43 




NGC 4636 


E 


-21 


.86 


203 


n 


2 


.21 


16. 


,76 




NGC 4648 


E 


-20 


.24 


220 


\ 





.83 


13. 


,34 




NGC 4649 


E 


-22, 


.51 


336 


n 


2, 


.34 


16. 


,77 




NGC 4660 


E 


-20, 


.13 


188 


\ 


0, 


.54 


12. 


,53 




NGC 4696 


BCG 


-24, 


.33 


254 


n 


2 


.44 


17. 


,77 


A3526-M1 


NGC 4697 


E 


-21 


.49 


174 


\ 





.41 


14. 


,13 




NGC 4709 


E 


-22, 


.32 


242 


n 


2 


.02 


16. 


,91 




NGC 4742 


E 


-19, 


.90 


109 


\ 





.51 


12. 


,43 




NGC 4874 


E 


-23, 


.49 


278 


n 


2, 


.99 


18. 


,98 




NGC 4889 


BCG 


-23, 


.73 


401 


n 


2, 


.84 


17. 


,80 


A1656-M1 


NGC 5017 


E 


-20 


.67 


184 


\ 





.99 


13. 


,30 




NGC 5061 


E 


-22 


.01 


186 


n 


1 


.39 


14. 


,06 




NGC 5077 


E 


-22, 


.07 


256 


n 


1, 


.96 


16. 


,07 




NGC 5198 


E 


-21, 


.23 


196 


n 


1, 


.33 


15. 


,19 




NGC 5308 


SO- 


-21, 


.26 


211 


\ 


0, 


.90 


13. 


,15 




NGC 5370 


SO 


-20, 


.60 


133 


\ 


1, 


.04 


15. 


,34 




NGC 5419 


E 


-23 


.37 


333 


n 


2, 


.65 


17. 


,53 




NGC 5557 


E 


-22 


.62 


254 


n 


1 


.82 


15. 


,58 




NGC 5576 


E 


-21, 


.31 


183 


n 


1 


.21 


14. 


,38 




NGC 5796 


E 


-21, 


.98 


288 


A 


1, 


.02 


14. 


,40 




NGC 5812 


E 


-21, 


.39 


200 


\ 


0, 


.84 


14. 


,27 




NGC 5813 


E 


-22, 


.01 


239 


n 


1, 


.89 


16. 


,32 




NGC 5831 


E 


-21, 


.00 


164 


\ 


0, 


.85 


14. 


,41 
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Table 1 — Continued 









a 




7 






Galaxy 


Morph 


My 


(km/s) 


P 


loefoc) 


(V-Band) 


Alt-ID 


NGC 5838 


SO- 


-20.51 


266 


\ 


1.03 


13.61 




NGC 5845 


E 


-19.98 


234 


\ 


1.14 


13.81 




NGC 5898 


E 


-21.65 


218 


A 


1.43 


15.41 




NGC 5903 


E 


-21.90 


198 


n 


2.17 


17.07 




NGC 5982 


E 


-21.97 


240 


n 


1.80 


15.62 




NGC 6086 


BCG 


-23.11 


336 


n 


2.53 


17.26 


A2162-M1 


NGC 6166 


BCG 


-23.80 


310 


n 


3.17 


19.32 


A2199-M1 


NGC 6173 


BCG 


-23.59 


278 


n 


2.32 


16.72 


A2197-M1 


NGC 6278 


SO 


-20.81 


150 


\ 


0.99 


13.97 




NGC 6340 


so 


-19.46 


144 


\ 


0.91 


14.54 




NGC 6849 


so 


-22.78 


216 


n 


1.98 


16.81 




NGC 6876 


E 


-23.58 


234 


n 


2.17 


17.02 




NGC 7014 


BCG 


-22.18 


263 


n 


1.83 


15.54 




NGC 7052 


E 


-22.35 


271 


n 


2.29 


16.19 




NGC 7213 


Sa 


-21.71 


182 


n 


1.83 


15.88 




NGC 7332 


SO 


-19.62 


124 


\ 


0.67 


12.78 




NGC 7457 


SO- 


-18.62 


69 


\ 


0.43 


15.86 




NGC 7578B 


BCG 


-23.41 


214 


n 


2.06 


16.19 


A2572-M1 


NGC 7619 


E 


-22.94 


322 


n 


2.03 


15.90 




NGC 7626 


E 


-22.87 


276 


A 


1.66 


14.98 




NGC 7647 


BCG 


-23.97 


282 


n 


2.28 


17.14 


A2589-M1 


NGC 7727 


Sa 


-21.19 


196 


A 


0.48 


14.11 




NGC 7743 


S0+ 


-20.18 


84 


\ 


1.03 


14.07 




NGC 7785 


E 


-22.08 


245 


n 


1.32 


15.28 




IC 0115 


BCG 


-22.67 




n 


2.45 


17.24 


A0195-M1 


IC 0613 


BCG 


-22.27 


262 


n 


2.05 


16.25 


A1016-M1 


IC 0664 


BCG 


-22.86 


336 


n 


2.07 


15.81 


A1142-M1 


IC 0712 


BCG 


-23.29 


345 


n 


2.69 


17.68 


A1314-M1 


IC 0875 


SO 


-20.21 




\ 


1.01 


13.45 




IC 1459 


E 


-22.51 


306 


n 


1.94 


15.39 




IC 1565 


BCG 


-22.99 


303 


n 


1.65 


16.86 


A0076-M1 


IC 1633 


BCG 


-23.91 


355 


n 


2.43 


16.66 


A2877-M1 
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a 










Galaxy 


Morph 


My 


(km/s) 


P 


losi(vc) 


(V-Band) 


Alt-ID 


IC 1695 


BCG 


-23.90 


364 


n 


2.36 


16.68 


A0193-M1 


IC 1733 


BCG 


-23.43 


301 


n 


2.68 


17.63 


A0260-M1 


IC 2738 


BCG 


-22.18 


275 


\ 


1.57 


16.15 


A1228-M1 


IC 4329 


BCG 


-23.95 


275 


n 


2.34 


16.26 


A3574-M1 


IC 4931 


BCG 


-23.47 


273 


n 


2.42 


16.86 


A3656-M1 


IC 5353 


BCG 


-22.64 


262 


n 


2.04 


16.37 


A4038-M1 


UGC 4551 


SO 


-19.78 


167 


\ 


0.82 


15.00 




UGC 4587 


so 


-20.77 




\ 


1.05 


14.79 




UGC 6062 


so 


-20.34 




\ 


1.01 


14.66 




VCC 1199 


E 


-20.34 




\ 


0.54 


15.42 




VCC 1440 


E 


-20.34 




\ 


0.54 


15.62 




VCC 1545 


E 


-20.34 




\ 


0.54 


17.49 




VCC 1627 


E 


-20.34 




\ 


0.54 


15.68 




ESQ 378-20 


SO 


-20.97 




\ 


1.06 


13.11 




ESQ 443-39 


so 


-20.93 




\ 


1.06 


14.97 




ESQ 447-30 


so 


-21.17 




\ 


1.00 


13.99 




ESQ 462-15 


E 


-22.83 


305 


\ 


1.19 


14.75 




ESQ 507-27 


SO 


-20.89 


203 


\ 


1.08 


14.25 




ESQ 507-45 


SO 


-23.28 


311 


A 


1.81 


14.95 




MCG 11-14-25A 


E 


-19.08 




n 


1.38 


15.77 




MCG 08-27-18 


SO 


-20.03 


89 


\ 


1.07 


15.09 




A0119-M1 


BCG 


-24.01 


294 


n 


2.81 


18.52 




A0168-M1 


BCG 


-23.12 


345 


n 


2.00 


16.90 




A0189-M1 


BCG 


-21.89 


230 


\ 


1.51 


17.56 




A0261-M1 


BCG 


-22.95 




\ 


1.64 


18.22 




A0295-M1 


BCG 


-23.11 




n 


2.63 


17.81 




A0376-M1 


BCG 


-23.60 


276 


n 


2.64 


18.10 




A0397-M1 


BCG 


-23.42 


289 


n 


2.70 


17.80 




A0419-M1 


BCG 


-21.79 




\ 


1.58 


17.08 




A0496-M1 


BCG 


-24.28 


273 


n 


2.61 


18.14 




A0533-M1 


BCG 


-22.68 




n 


2.28 


17.08 




A0548-M1 


BCG 


-22.75 


220 


n 


2.22 


17.12 
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Galaxy 


Morph 


Mv 


a 

(km/s) 


P 


1 

loff(pc) 


(V-Band) Alt-ID 


A0634-M1 


BCG 


-22.70 


245 


n 


2.17 


17.16 


A0912-M1 


BCG 


-22.24 




n 


1.63 


16.36 


A0999-M1 


BCG 


-22.45 


272 


n 


2.29 


17.03 


A1020-M1 


BCG 


-22.65 


345 


n 


2.32 


16.87 


A1631-M1 


BCG 


-23.34 


249 


n 


2.12 


16.49 


A1831-M1 


BCG 


-23.51 




n 


2.84 


18.58 


A1983-M1 


BCG 


-22.35 


270 


A 


1.73 


15.50 


A2040-M1 


BCG 


-23.46 


223 


n 


2.28 


17.41 


A2052-M1 


BCG 


-23.04 


216 


n 


2.46 


18.53 


A2147-M1 


BCG 


-23.16 


278 


n 


2.90 


19.03 


A2247-M1 


BCG 


-22.66 


209 


\ 


1.60 


20.06 


A3144-M1 


BCG 


-22.28 




n 


2.28 


16.79 


A3376-M1 


BCG 


-23.29 


300 


n 


3.11 


18.93 


A3395-M1 


BCG 


-24.23 


276 


n 


2.52 


18.09 


A3528-M1 


BCG 


-24.30 


434 


n 


2.61 


17.54 


A3532-M1 


BCG 


-24.58 




n 


2.51 


17.49 


A3554-M1 


BCG 


-23.99 




n 


2.62 


18.36 


A3556-M1 


BCG 


-23.65 




n 


2.48 


16.95 


A3558-M1 


BCG 


-24.92 


275 


n 


3.12 


19.29 


A3562-M1 


BCG 


-24.32 


236 


n 


2.84 


18.87 


A3564-M1 


BCG 


-22.68 




n 


2.12 


16.63 


A3570-M1 


BCG 


-22.54 


268 


n 


2.01 


16.19 


A3571-M1 


BCG 


-25.30 


325 


n 


3.03 


19.39 


A3677-M1 


BCG 


-22.21 




n 


2.14 


16.49 


A3716-M1 


BCG 


-23.75 


247 


n 


2.56 


17.99 


A3736-M1 


BCG 


-23.98 




n 


2.70 


17.93 


A3747-M1 


BCG 


-22.65 


232 


n 


2.00 


15.97 



Note . — Derivation of the parameters listed are presented in iLauer et al. 
(|2007al ). Bulge luminosities are given for SO and spiral galaxies. Velocity 
dispersion are provided by the "Hyperleda" database. The profile type, P, is 
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\ = power-law, A = intermediate form, and n = core. BCG identificat ions 
are referred to their hosting Abell clusters; see lPostman &: Laueii (|l995l ) for 
details. 
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Table 2. Core Galaxies with Measured Black Hole Masses 



Galaxy 


My 


log(r^/pc) 


I, (V) 


log(M./M0) 


M, Source 


NGC 1399 


-22.07 


2.23 


16.76 


8.95 


1, 2 


NGC 3379 


-21.14 


1.72 


15.59 


8.00 


3 


NGC 3608 


-21.12 


1.31 


15.05 


8.28 


4 


NGC 4261 


-22.26 


2.31 


16.09 


8.72 


5 


NGC 4291 


-20.64 


1.63 


15.29 


8.49 


4 


NGC 4374 


-22.28 


2.11 


15.67 


9.00 


6, 7 


NGC 4473 


-21.16 


1.73 


15.40 


8.04 


4 


NGC 4486 


-22.71 


2.65 


17.25 


9.48 


8 


NGC 4649 


-22.51 


2.34 


16.77 


9.30 


4 


NGC 7052 


-22.35 


2.29 


16.19 


8.52 


9 


IC 1459 


-22.51 


1.94 


15.39 


9.18 


10, 11 



Note. — Black hole mass references are 1) Houehton et al. 


(2006^. 2) Gebhardt et al. (20071. 3) Gebhardt et al. ( 


2000b 


), 4) 


Gebhardt et al.l (20031. 5) Ferrarese et al 


. (1996). 6) ] 


Bower et al. 


(19981. 7) Macieiewski & Binnev 


J2001), 8) 


Macchetto et a] 


. (1997). 9) 


van der Marel & van den Bosch ( 


19981). 10) ' 


V^erdoes Kleiin et al. ( 


2000), 


and 11) Cappellari et al. (2002). 


For galaxies with two references, the 



black hole mass is an average value. 
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Fig. 1. 

function of My. Galaxies are drawn from the 



Black hole masses for all black ho les with direct niass d eterminations are plotted as a 

I2OO2I ) sample with augmentations as 



Tremaine et al 



described in the text. The solid line is the lHaring &: Rixl (j2004l ) relationship (equation @) between 
M, and galaxy mass transformed to luminosity using M/Ly ~ Ly"^^ with zeropoint M/Ly = 6 at 
My = —22. A symmetrical least-squares fit to all data points is shown as the dashed line (equation 
[2]), and a fit to just the galaxies with My < — 19 is shown as the dotted line (equation [3]) . 
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Fig. 2. — M,{L) versus M,((t) the sample galaxies that have a measurements. Power-law galaxies 
are plotted as red dots, core galaxies are blue dots, "intermediate" galaxies are plotted as small 
open circles, and BCGs with resolved cores are plotted as green squares. Galaxies with large circles 
have directly determined black hole masses; however, the predicted rather than observed M, values 
are still plotted. The M, — a relationship is that of Tremaine et al. (2 002). The asy mmetric error 
bar in the horizontal direction shows the change in predicted M, if the (|Wvithdl2006l ) log-quadratic 
M, — cr relationship is used instead. M,(L) is the average of the minimum and maximum predictions 
for a given L from the three M, — L relationships in Figure [H with the error bars showing the range 
of the predictions. 
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Fig. 3. — The relationship between central velocity dispersion, cj, and L for the sample is plotted. A 
fit to just the core galaxies and BCGs (solid line; equation[5]) gives L ~ cj^, a much steep relationship 
then the standard L ~ Faber-Jackson relationship, and L ^ a'^ for the power-law galaxies alone 
(dashed line; equation [6]) . It is this change in slope that leads to conflicting predictions for M, 
from the M, — L and M, — a relation for the most luminous galaxies. Core galaxies with directly 
measured black hole masses are circled. 
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100 

o (km/s) 



Fig. 4. — Cusp radius, r^, is plotted as a function of stellar velocity dispersion. The power-law 
galaxies are now plotted as triangles to indicate that their cusp radii are only upper limits. The 
solid line is the fitted relationship (equation [9]) between and a for core galaxies. The figure shows 
that is a steep function of a. If M, ~ cr'^ as equation ([T]), the observed empirical relationship 
between and ex (solid line) implies that ~ 
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Fig. 5. — Cusp radius, r^, is plotted as a function of total galaxy luminosity. Power-law galaxies are 
plotted as V to indicate that their values are as upper limits. The solid line shows the best-fit 
relationship between and My for core galaxies (equation [TO]) . The figure shows that varies 
nearly linearly with L. If M, ~ Ly^ as in equation ([H), the observed empirical relationship between 
and L (solid line) implies that ~ JVf0.96±o.09_ 

Note that for core galaxies the range in at any given L is smaller than it is at a given a. 
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(parsec) 

Fig. 6. — Cusp brightness, /x^, is plotted as a function of cusp radius, r-y. The power-law triangles 
are rotated and shown as arrows to reflect that the points are only upper limits for both Xy and r-y. 
The tight relationship between and (equation I15|) means that either can serve for the other 
in the context of relating core structure to M,, L or a. 



- 52 - 




Fig. 7. — Black hole mass versus core size, r^, for the 11 core galaxies that have measurements. 
The red line is the symmetric fit between M, and provided by equation (j20l) . while the blue line 
gives the fit presented in equation (j2ip . which assumes that is the independent variable. 
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Fig. 8. — The four panels plot M, predicted from as function of a and My for core galaxies, 
where green symbols are BCGs, and blue symbols are the remaining core galaxies. The red lines 
give M, values predicted either from the M, — a relationship (equation [T]) or the M, — L relationship 
(equation [4]) . The upper panels give M, predicted by through equation ([20l) . which was derived 
by a symmetrical fit to the points in Figure [71 The bottom panels, in contrast, use equation (I2ip . 
which was derived assuming as the independent variable. Both equations typically predict M, 
in excess of the predictions from the M, — a relationship. 




108 



Mbh(^) (Ms) 



109 
Mbh(L) (MJ 



101 



Fig. 9. — The two panels plot core mass as function of M,((t) and M,{L) for core galaxies, where 
green symbols are BCGs, and blue symbols are the remaining core galaxies. The red lines give the 
mean — M, relationships inferred by combining either the M, — a relationship (equation [1]) with 
the — a (equation [T7|) , or the M, — L relationship (equation (H) with the My — L relationship 
(equation [T6|) . 
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Fig. 10. — Effective radius as a function of luminosity for the galaxy sample. The steepening of the 
Rg — L relationship sets in at My < —22, similar to the luminosity at which the velocity dispersion 
starts to "plateau" in Figure [3l We attribute this to a progressive change in the character of "dry 
mergers" at higher galaxy masses. The shallow line is defined by a fit to power-law galaxies only 
(equation I25p . while the steep line is a fit to core galaxies with My < — 21 (equation [26]) . Power-law 
galaxies and core galaxies, however, have similar at My ~ —21 where the transition between 
the two forms takes place. 
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Fig. 11. — The log cumulative density of black holes above a given mass versus logM, for different 
ma ss functions. The red c urve in both panels is derived from a lightbulb model for quasars applied to 
the lRichards et al.l (120051 ) best-fit quasar luminosity function from SDSS, evaluated at z = 2.5. The 
model assumes they radiate at their Eddington Luminosity with a duty fraction / = . 03 (se e text). 
The pink curve in both panels is the BH mass function produced by the lHopkins et al.l (|2006l ) model 
for quasar, described in detail in the text. The cya n curve i n both panels is t he BH mass function 
obtained by augmenting th e iBlanton et al.l (120031) b est-fit ISchechteii (|l976l ) luminosity function 



for SPSS galaxies with the iPostman &: Laueii (jl995l ) brightest cluster galaxies and calibrated by 
the iHaring &: Rixl ([200J) relation. The solid and dotted lines above the curve sho w the effect of 
cosmi c scatter of 0.25 and 0.50 (respectively) about the mean relation (derived by iHaring &: Rix 
(|2004l )) between BH mass and galaxy mass. The green curve (in the left panel ) is a BH mass 



fun ction predicted f rom th e SDSS velocity dispersion function (jSheth et al 



20031) augmen t ed by 



the iBernardi et al.l ()2006al ) high-dispersion sample and calibrated by the iTremaine et al.l ()2002l ) 
M,{a) predictor for zero cosmic scatter. The solid and dotted lines above the curve show the effect 
of a cosmic scatter of 0.15 and 0.30 in the decimal log of the BH mass about the mean Tremaine 
relation. The dark blu e curve in the right panel illustrates the same dispersion function calibrated 
instead by the IWvithd (j2006l ) M,{a) predictor. 
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Fig. 12. — SPSS r "mo del magnitudes" are compared to two separate Rc magnitudes for the 



Postman &: Laueii (|l995l ) BCGs inc common, where the SDSS values are transformed assum- 



ing Rc = i ~ 0-3 1. Red points show the difference between SDSS model r magnitudes and the 
Laine et ah ( 2002) Rg total BCG l uminosity versus effective radius, Re, which is derived from r^/^ 



fits to the lPostman &: Laueii (119951) surface brightness profiles. The blue points are the same exer- 
cise, but with the lPostman &: Laueii (|l995l ) maximum-aperture magnitudes used instead. Clearly, 
the larger a BCG is, the more the present total luminosity disagrees with the SDSS value. The 
maximum- aperture magnitudes are not intended to be interpreted as a total magnitude, but pro- 
vide a model-independent lower limit on its value. Since the maximum-aperture magnitudes are 
brighter than the SDSS magnitudes for most BCGs, this clearly shows that the SDSS values cannot 
be regarded as total luminosities. 
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Fig. 13. — Surface brightness profiles from iPostman &: Laueii (|l995l ) (black) are compared to 
"profmean" SDSS profiles (red) for three BCGs. The SDSS r band magnitude s are t ransformed 
assuming Rq = r — 0.31. The SDSS profiles all fall below the lPostman &: Laueii (|l995l ) profiles at 
large radii, consistent with excessive sky subtraction by the SDSS pipeline. 
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Fig . 14. — Histograms of es timated total My are shown for three BCG samples. Magnitudes for 
the Postman &: Lauer volume- l i mited BCG sample (blue) are based on r^/^ fits to -Rc-band 

surface photometry. The lMiller et al.l (|2005l ) BCG sample i s based on SDSS model r magnitudes, 
and has typical luminositie s one magnitude sinaller t han the lPostman &: Laueii (|l995l ) sample. The 
histogra m of the sub s et of Postman k. Laueii (jl995l ) BCGs observed by SDSS (green) agrees well 



with the lMiller et al.l (|2005l ) sample when SDSS r model magnitudes are used instead to estimate 
total My, yet we argue that these magnitudes are strongly affected by excessive sky subtraction. 
This concordance implies that the C4 BCGs, are also likely to have had their total luminosities 
under-estimated . 
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Fig. 15.— The archived 2MASS J b and image of NGC 2832 the BCG in A0779, is compared to a 
portion of the R band image used bv lPostman &: Laueii (|l995l ) to derive surface photometry profile 
shown in the next figure. The stretch has been set to be the same for both images. The R band 
image has been binned to a 0"91 pixel scale to roughly match the I'.'O scale of the J band image. 
The sky level of the 2MASS image is effectively 26 x brighter in J, taking the observed R — J color 
of the galaxy into account. The J band image is clearly considerably shallower than the R band 
image, and the envelope of the galaxy disappears into the noise at radii where it is still clearly 
present in the R band. 
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Fig. 16. — Major axis surface J photometry profi les (red) derived f r om 2M ASS J band archive 
images are compared to the R profiles derived by iPostman &: Laueii (|l995l ) (black) for the three 
BCGs shown in Figure [T3l The ~ 0.4% error in the 2MASS sky levels gives the large error bars. 
The last J band isophotes fall ~ 7 magnitudes below the sky, and thus are less significant than 
the errors in the sky levels. The R and J profiles agree within the errors, and total luminosities 
estimated by r^/^-law fits to the even the J band profiles are considerably larger than the 2MASS 
XSC apparent luminosities. 
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Fig. 17. — Residuals about the mean r;, — My relationship are plotted for core galaxies as a function 
of the logarithmic envelope slope /?. At any given My, excessively large cores (positive residuals) 
correspond to higher (5 (steeper envelopes) while excessively small cores correspond to small (5. 
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Fig. 18. — The — L and rf, — L relationships for core galaxies are compared. The hnes are the 
mean relationships for My < —21. The — L relationship has smaller scatter for core galaxies 
with My < —21. Power-law galaxies are plotted the same in both panels, with upper limits on 
to be the same as upper limits on r^. 



